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Point-contact Andreev reflection spectroscopy (PCAR) has proven to be one of the most powerful 
tools in the investigation of superconductors, where it provides information on the order parameter 
(OP), a fundamental property of the superconducting state. In the past 20 years, successive im- 
provements of the models used to analyze the spectra have continuously extended its capabilities, 
making it suited to study new superconductors with "exotic" properties such as anisotropic, nodal 
and multiple OPs. In Fe-based superconductors, the complex compound- and doping-dependent 
Fermi surface and the predicted sensitivity of the OP to fine structural details present unprecedent 
challenges for this technique. Nevertheless, we show here that PCAR measurements in Fe-based 
superconductors carried out so far have already greatly contributed to our understanding of these 
materials, despite some apparent inconsistencies that can be overcome if a homogeneous treatment 
of the data is used. We also demonstrate that, if properly extended theoretical models for Andreev 
reflection are used, directional PCAR spectroscopy can provide detailed information not only on 
the amplitude and symmetry of the OPs, but also on the nature of the pairing boson, and even give 
some hints about the shape of the Fermi surface. 

PACS numbers: 74.50.-fr , 74.70. Dd, 74.45. +c 



I. INTRODUCTION 

The order parameter (OP in the following) is the most 
important quantity in a superconductor, since it deter- 
mines how the charge carriers couple to form Cooper 
pairs. In the standard BCS theory for superconductivity 
[J , the pairing potential is always attractive, does not de- 
pend on energy, is isotropic in space and is due to the in- 
teraction between electrons and phonons (in other words, 
phonons are the "mediating boson"). The amplitude of 
the OP, A, is also the "energy gap" in the superconduc- 
tor, i.e. the distance, in energy, between the condensate 
of Cooper pairs and the first single-particle excitation. A 
dependence on energy of the OP must be added, as in the 
so-called Eliashberg theory [5] , to properly describe those 
superconductors, like Pb, where the interaction between 
electrons and phonons is particularly strong. In the past 
decades, many superconductors have been discovered in 
which the OP is anisotropic (i.e. its amplitude shows 
some angular dependence) or can even change sign to 
become negative in some directions (as in the so-called 
d-wave symmetry, common in copper-oxide superconduc- 
tors); in some cases, like the well-known MgB2, multiple 
OPs exist, related to different sets of bands that form 
separated "sheets" of the Fermi surface. 

The recently discovered Fe-based superconductors 
present most of these unusual properties at the same 
time. The coupling of electrons with the mediating bo- 
son (whatever its nature is) seems to be intermediate to 
strong; the Fermi surface (FS) is made up of separated 
sheets, which are mostly 2D-like; there are probably mul- 
tiple OPs that, in most cases, are expected to be isotropic 
but with opposite sign on hole-like and electron-like FS 
sheets [3]. In some conditions, OPs that display a sign 



change on the same sheet of the FS are predicted to be- 
come energetically favorable [H |3] . 

As usual when a new class of superconductors is dis- 
covered, the experimental study of the OP in Fe-based 
superconductors has immediately attracted the attention 
and expectations of the scientific community. Determin- 
ing the number, the amplitude and the symmetry of the 
OPs is a necessary step to assess the coupling mechanism 
and to understand the details of how superconductivity 
occurs in these materials. Point-contact spectroscopy is 
one of the techniques that have been more useful in the 
investigation of the OP in these materials. Despite its 
apparent simplicity (it just consists in measuring the dif- 
ferential conductance of a very small, in principle point- 
like contact between a normal metal and the supercon- 
ductor under study) it is a powerful and versatile, yet 
inexpensive, research tool generally based on a combi- 
nation of quantum effects: the tunnel effect (dominant 
in high-resistance contacts) and the Andreev reflection 
(dominant in small-resistance contacts). Here we will 
concentrate on the latter. To explain in a simple way 
what Andreev reflection is, let's imagine an incident elec- 
tron with total energy E to travel towards the interface 
in the N (normal) side of a ideal, barrierless S — N junc- 
tion {S = superconductor). If it finds vacant electronic 
states at the same energy on the S side, the electron will 
simply propagate in S. If instead its energy is smaller 
than the energy gap A in S*, it cannot propagate in 5, so 
it forms a Cooper pair with another electron excited from 
the Fermi sea in S. This leaves a hole in S that crosses 
the interface and travels backwards in N (see [SHE] for 
more accurate descriptions). From the point of view of 
N , things happen as if the electron was reflected as a 
hole. This "reflection" causes a doubling of the conduc- 
tance with respect to a N — N junction in which Andreev 
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reflection cannot take place. This doubling occurs only 
as long as the energy of the injected electron (equal to 
eV if V is the voltage difference between N and S) is 
smaller or equal to A. Thus, a measurement of the junc- 
tion conductance as a function of V allows determining 
the gap, i.e. the amplitude of the OP. In unconventional 
superconductors, the shape of the conductance curves 
can provide information not only on the amplitude of the 
OP, but also on its angular dependence (i.e. symmetry 
in reciprocal space) and, in particular conditions, on the 
shape of the Fermi surface and on the spectral function 
of the boson that mediates the superconducting pairing, 
as we will discuss in section Hi CI 

Point-contact Andreev-reflection (PCAR) spec- 
troscopy measurements in Fe-based superconductors 
have been carried out by various groups. Early mea- 
surements suffered from the unsatisfactory quality of 
the first samples and gave indeed contradicting and 
sometimes unreproducible results. This might have 
given the impression that PCAR does provide reliable 
results in these compounds. We will show in section |III| 
that, instead, when PCAR experiments in good-quality 
samples are carried out with all the necessary cautions 
by different groups and even in samples of different 
forms (single crystals, polycrystals, films) they do 
provide surprisingly consistent results. The sometimes 
contradicting conclusions arise from the interpretation 
of the data rather than from the data themselves. We 
will thus show that, in the most studied 1111 and 122 
compounds, the available PCAR results agree rather 
well on the amplitude of the OPs, on their temperature 
dependence and on their symmetry (the s± one, with 
isotropic OPs with sign reversal between holelike and 
electronlike FS sheets being the most likely, at least 
at optimal doping). In the same section we will also 
present and discuss recent PCAR results in the 122 and 
22426 families that seem to indicate a nodal, or strongly 
anisotropic OP. If confirmed, these findings can serve 
as a test for the theories that predict the occurrence of 
anisotropic and nodal symmetries in suitable conditions, 
related to fine details of the lattice structure. 

The problem of extracting the electron-boson spec- 
trum from PCAR measurements in these materials will 



be dealt with in section IV We will show that, with the 



combined use of Eliashberg theory and of the models for 
PCAR conductance, it is possible to extract the charac- 
teristic energy of the mediating boson and also to pose 
some constraints on the shape of the relevant spectrum. 
The strong electron-boson coupling also explains some 
anomalies in the PCAR conductance curves, namely the 
excess conductance at high energy, that greatly compli- 
cates the normalization of the conductance curves and 
sometimes prevents a good fit of the curves with models 
based on constant, BCS-like OPs. 

Finally, in sectionjV] we will show that the values of the 
gap ratios 2Ai/kBTc collected from PCAR measurements 
in various Fe-based superconductors seem to follow a uni- 
versal trend as a function of the critical temperature (T^) 



of the samples. Both 2Ai/kBTc and 2A2/ ksTc remain 
approximately constant down to ~ 30 K, then tend 
to increase in compounds with lower - even though 
further measurements in the low- Tc region would be nec- 
essary to confirm this trend. We will propose a possible 
explanation of this behavior that essentially relies on the 
fact that the mediating boson is an electronic excitation 
(for example, spin fiuctuations) and thus is subject to a 
"feedback effect" of the condensate on the electron-boson 
spectral function. 



II. MULTIBAND 3D MODELS FOR 
DIRECTIONAL ANDREEV-REFLECTION 
SPECTROSCOPY 

A. The "standard" ID and 2D BTK models 

We already presented in a recent review [9] the ex- 
plicit expressions for the normalized conductance of a 
point contact between a normal metal (N) and a super- 
conductor (S) in the case of a 2D geometry of the Fermi 
surfaces of the latter. Let us briefly summarize these 
equations here, since in subsection |II C| we will use them 
as the starting point to extend the Blonder-Tinkham- 
Klapwijk (BTK) model to particular 3D shapes of the 
FS and to particular directions of current injection. In 
the following we will call "2D BTK model" the gener- 
alization of the original, "ID" BTK theory [3], intro- 
duced in 1996 by Y. Tanaka and S. Kashiwaya [7]. In 
their article they first wrote the normalized conductance 
G{E) = {dI/dV)T~^s/{dI/dV)m at T = (where Ins 
and /nn are the currents flowing through the interface 
when the material under study is in the superconduct- 
ing and in the normal state, respectively) as a func- 
tion of the two quantities Nq{E) = E/^/E"^ — A"^ and 
Np{E) = A/VE^ - A2, whose real parts are the BCS 
quasiparticle and pair density of states, respectively. 

Here we will immediately generalize these expressions 
to the strong-coupling regime, where the order param- 
eter A is a complex function of energy resulting from 
the solution of the Eliashberg equations, and to the 
presence of inelastic quasiparticle scattering processes 
at the N/S interface. In the latter case it has been 
shown [ini that the resulting (and dominant) extrinsic 
reduction of the quasiparticle lifetime can be properly 
taken into account by including in the model an imagi- 
nary part of the energy described by the broadening pa- 
rameter r, i.e. E —i- E + iV. As a consequence, from 



iT) 



- A{EY and 
Note that we 



now on: A^q(£:,r) = (S -f 

7Vp(i?,r) = A{E)/^{E^iTY-A{EY 
have allowed A to depend on energy for the sake of gener- 
ality, even though wc will use this dependence explicitly 



only in section IIP 



The main steps in the calculation of the normalized 
conductance G{E) consist in the determination of the 
barrier transparency when the material under study is 
in the normal state (tn) and in the superconducting 



3 



state {(Js{E)). These quantities are actually transmission 
probability distributions and are obtained by solving the 
Bogoliubov-De Gennes equations with the only condition 
that the component of the wavevector k parallel to the 
interface is conserved in all processes, for any direction 
of the incoming electron from the N side. In the pure 
2D case here discussed, where we suppose the interface 
to coincide with a plane parallel to the z axis, the incom- 
ing electrons lie in the xy plane and have a wave vector 
kp.N which forms an angle 0n with respect to the uni- 
tary vector n normal to the interface. In principle, if the 
Fermi velocity is different in the N and S side, the charge 
carriers experience a "refraction" when crossing the in- 
terface (i.e. the wavevector of the quasiparticles in the S 
side of the junction kp s could be different from kp n) so 
that the direction of propagation of quasiparticles in the 
S side could form an angle 9s different from 9^ with re- 
spect to the normal to the interface. In these conditions, 
electrons incident on the interface with suitable angles 
can experience a normal reflection even in absence of a 
barrier potential Uq located at the interface, and the nor- 
mal transmission coefficient tn is no longer identically 1. 
In the most general case when Uq is present, the normal 
transparency of the barrier is: 

4(kp^N • n)(kp,s • n) 



(kp^N • n + kp3 • n)2 -f- AZ^k 
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where Z — Uq/Hvp is the dimensionless parameter which 
accounts for the potential barrier. 

From now on we will restrict the analysis to the case 
kp N=kp g = kp from which it follows 9j^=9s=9. As a 
matter of fact, while properly taking into account the 
Fermi velocity mismatch at the interface in the 2D case 
is a non-straightforward but feasible task (see [9]), it be- 
comes an extremely complex problem in the full 3D case 
we will discuss in subsection |II C[ when the true shape of 
the FS will be considered. Under the previous restrictive 
hypothesis the normal transparency becomes: 



MO) = 



cos 



'9 + Z^' 



(2) 



It is clear that, for Z = 0, the normal transmission 
probability is identically 1 for any direction of the in- 
coming electron, i.e. for — 7r/2 < 9 < Tr/2. When Z ^ 0, 
instead, the barrier transparency depends on the direc- 
tion of the incoming electron in the N side; for Z = 10 
(tunneling regime) the transmission probability is always 
small and highly directional. It can be shown that the 
relative transparency of the barrier in the superconduct- 
ing state is given by [7]: 



1 + M9ME)\^ + {t^{9)-IM E) 
\l + [M9)-l]j{E)^\^ 



2|2 



(3) 



where -/{E) = [Nq{E) - 1] /Np{E) and finally the nor- 
malized conductance at T = is: 



G2T>{E) — 



JZi as{E,9)M9)cos^ 

2 

/+! rN(6') cos 9d0 



(4) 
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FIG. 1: (color online) ; (a) Theoretical PCAR normalized con- 
ductance of a single-band s-wave superconductor as function 
of reduced energy _E/A at T = and F = for different Z 
values from the pure Andreev-reflection regime {Z = 0) to 
the tunneling one {Z = 10). The solid curves are calculated 
by using the ID BTK model, while the dashed ones by using 
the 2D BTK one (see text), (b) Temperature dependency of 
the theoretical PCAR normalized conductance of panel (a) for 
A = 3 meV, F = and Z — 0.3. Also in this case solid curves 
represent the results of ID BTK model and dashed ones the 
results of 2D BTK one. In the inset the BCS temperature 
dependency of the gap is shown. The curves in (a) and (b) 
are vertically offset for clarity. 



Both at numerator and denominator of equation [4] the 
integration extremes are — 7r/2 and tt/2 since the current 
injection occurs in the whole positive half-space for which 
< cos 6* < 1. The standard BTK theory [H is the ID 
limit of this model, in which all the electrons from the N 
side are injected in the same direction perpendicular to 
the interface. Obviously, the normal-state transmission 
probability in this case is simply obtained by putting 
9 — in equation [2] i.e. tn = 1/(1-1- Z'^). By doing the 
same also in equations [3] and |4] the expression for the ID 
BTK conductance Giu{E) — asiE,0) is thus recovered. 
Instead, the 2D BTK model itself is the limit of the more 



general 3D model (see section II C I if the gap is isotropic 



and the FS is spherical (i.e. the system has a rotational 
symmetry around the axis normal to the interface). 

In Figure [h a) a comparison of the normalized conduc- 
tances Giu(E) (sohd hues) and G2d{E) (dashed lines) 
is shown at T = and F = as a function of the nor- 
malized energy _B/A for different values of the barrier 
parameter Z. Looking at the form of t^{9) (equation [2|, 
one expects the ID BTK normalized conductance to co- 
incide with the 2D one in two opposite conditions: when 
Z = Q (i.e. when tn = 1 and gs{E) ^ 1 + \i{E)\'^) 
and when Z — > 00 (i.e. when tn — > and cts{E) — >■ 
(1 - \-i{Ef\'^)/\l - ^{Ef\'^). This is clearly shown in 
the curves at Z = 10 and Z = 0. The maximum devia- 
tions between the predictions of the two models occur for 
Z ~ 0.3 — 0.5 which, incidentally, is the range of values 
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commonly observed in experiments in Fe-based supercon- 
ductors. It is very easy to demonstrate that, in practice, 
each G2d{E) curve almost coincides with the Gio{E) 
calculated for a properly enhanced Z value and, thus, 
the ID BTK model can be efficiently used instead of the 
2D one provided that the gap is isotropic, the FS spheri- 
cal and one is not interested in the precise determination 
of the Z values. 

All the equations seen up to now are calculated at 
T — 0. The effect of a finite temperature T can be prop- 
erly taken into account by simply convolving the normal- 
ized conductance with the Fermi function at that tem- 
perature, as shown in detail in [3] and by including in 
the calculations the proper temperature dependence of 
the gap. Figure [T|b) shows the temperature dependence 
of the normalized conductances Gid(-E) (solid lines) and 
G2d{E) (dashed lines) for A = 3.0 meV, F = 0, Z = 0.3 
and Tc = 20 K. The BCS temperature dependence of 
the gap A(r) used in the calculations is shown in the 
inset. This kind of curves is what is expected in the best 
PCAR experiments (F = 0) on a single-band isotropic 
BCS superconductor. 



B. Symmetry of the order parameter 

In order to understand the profound changes occurring 
in the normalized conductance of a point contact when 
the symmetry of the order parameter is not s wave, it is 
necessary to brieffy describe in greater detail the physi- 
cal processes that take place at the NS interface. When 
an electron coming from the N side reaches the potential 
barrier at the interface, four different processes can occur 
depending on the electron energy and on the height of the 
barrier: A) the electron is reflected as a hole in N (An- 
dreev reflection); B) the electron is specularly reflected 
as an electron in N (normal reflection); C) the electron 
is transmitted in S as an electron-like quasiparticle, ELQ 
(normal transmission); D) the electron is transmitted in 
S as a hole-like quasiparticle, HLQ (anomalous transmis- 
sion). Due to the charge conservation this last process 
results in two electrons "reflected" back in N for every 
incoming electron. The probabilities of these processes 
strongly depend on the energy of the incoming electron 
and on the Z parameter. In particular, processes C and 
D are forbidden for E < A (where there can be no trans- 
mitted quasiparticles) , process B is present at any energy 
when Z > and process A is also present at any energy 
(but is very small for E > A) only if Z values are not 
too large (in practice Z < 10). Equation |4] takes prop- 
erly into account all these probabilities, but only in the 
case of an isotropic order parameter, i.e. in s-wave sym- 
metry. When the OP is anisotropic, i.e. its amplitude 
and possibly its sign change as a function of the direc- 
tion in reciprocal space, the situation is more complex. 
As a matter of fact, the wavevector of the ELQ kp and 
that of the HLQ — kp have now different direction, yet 
still having the same component parallel to the interface. 



In other words, the two vectors in the S side form the 
angles 9 and —9 with respect to the unitary vector n, 
respectively [7]. In this situation and depending on the 
orientation of the anisotropic OP with respect to n the 
ELQs and the HLQs can feel different OPs. 

Let us extend the equations for the 2D BTK model to 
this case. Let us assume that the system has a transla- 
tional invariance along k^, i.e. the FS is a cylinder, the 
OP has a k dependence only in the k^ky plane (as it hap- 
pens in cuprates) and the current injection occurs in the 
same plane (a&-plane contact). We can now introduce the 
angle of misalignment a between the crystallographic a 
axis in S and the normal to the interface n. As a conse- 
quence of gap anisotropy and of this misalignment, ELQs 
feel the order parameter A{E, 0)+ = A{E, 9 — a), while 
HLQs feel A(i;,6')_ = A{E,-9-a). In this case the rel- 
ative transparency of the barrier in the superconducting 
state becomes [7]: 



as{E,9)^ 

l + T^i9)\j+{E,t 



(5) 

{M9)-l)\j+{E,9h.(E,9)\^ 



|1 + (tn(^) - lh+{E,9)j_{E,9)eMiVdW 



where 
and 



(E + 2T) - ^{E + »F)2 - \AiE:9)^ 
\AiE,9)±\ 



-i In 



AiE,9)+/\A{E,9) + \ 
A{E,9)^/\A{E,9U 



(6) 



(7) 



is the phase difference (function of 9) seen by the HLQs 
with respect to the ELQs. Of course, t^{9) remains the 
same (indeed it only describes processes B and C in a 
N-N junction) and thus, simply by introducing the new 
expression for as{E,9) (equation [s]) into equation [2] and 
performing the proper integration in 9 one can obtain the 
normalized conductance at T = for any symmetry of 
the OP in the fcx^y plane. This symmetry is of course 
expressed by the speciflc 9 dependence of the functions 
AiE,9)±. 

In the present review we are interested in the sym- 
metries that have been considered plausible for Fe-based 
superconductor. The s± symmetry, with isotropic OPs 
of different sign on the hole-like and on the electron-like 
FS sheets is the most likely for the majority of the 
systems. If interband interference effects [11] are negli- 
gible, this phase change is not detectable by PCAR and 
gives the same spectra as a multiband s-wave symme- 
try. In suitable conditions, that seem to be related to 
flne structural parameters [H [5], anisotropic or nodal 
symmetries are expected to become more energetically 
favorable. As a consequence, apart from the s-wave, 
we will discuss here the cases of fully anisotropic s- 
wave symmetry (with zeros in the gap) and d^^ -y2 -wave 
one. The two expressions for the OP are: Aan(£', ^)± = 
A(£;)cos4 [2{±9-a)] and Ad(-E,6')± = A{E) cos2{±9- 
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FIG. 2: (color online): (a) Reduced energy dependence of the normalized conductance of a single-band superconductor with 
full anisotropic s symmetry of the OP for different Z values at T = 1.5 K,a = 7r/4 and F = (upper panel) and at T = 1.5 
K,Q = 7I-/4 and F = A/4 (lower panel); (b) The same as in panel (a) but for a single-band superconductor with d^2_y2 symmetry 
of the OP; (c) The same as in panel (b) but for a = tt/S. The insets in all the panels represent the symmetry of the pair 
potential (red and blue colors mean positive and negative values, respectively) while the black arrows indicate the direction of 
current injection. 



a), respectively [2l[T2]- In the first case the pair potential 
is zero in four directions but never changes sign, while in 
the second case the sign change typical of d-wave sym- 
metry is present, as shovi^n in the insets of figure [2j 

Figure [2ja) reports the theoretical 2D AR normal- 
ized conductance curves obtained from the integration 
of equation |4] for fully anisotropic s-wave symmetry, at 
finite temperature (T — 1.5 K) and different Z values. 
Upper and lower panels refer to the cases F = and 
F = A/4, respectively. Figure |2][b) and (c) report the 
corresponding conductance curves obtained for a d-wave 
symmetry, when a = n/A and a = 7r/8, respectively. 
Without entering too much into details, the main mes- 
sage coming from these calculations is the following. For 
current injection in the ab plane, and when Z is suf- 
ficiently small, a maximum of the conductance at zero 
bias appears both in anisotropic s- and rf-wave symme- 
tries (top panels). These maxima however have a some- 
how different physical origin in the two symmetries. In 
fully anisotropic s symmetry, they come only from the 
angular integration over OPs that can assume any value 
between the maximum A.{E) and 0. In d-wave symme- 
try, the zero-bias maximum has the same origin only for 
Z ^ 0, while for finite Z values it considerably grows 
(becoming much higher than 2 in the tunneling regime) 
because of an additional effect, i.e. the constructive in- 
terference of ELQs and HLQs that feel a phase difference 
of the pair potential. This interference effect (described 



by the exp(z(prf) term in the denominator of equationjSJ) is 
the only responsible for the zero-bias conductance peak 
(ZBCP) present in the d-wave case for a > and large 
Z, a peak that is totally absent in the anisotropic s case. 

In principle, the different shape of the conductance 
curves for the two symmetries at low Z shown in the up- 
per panels of figure [2] would allow one to identify the true 
OP symmetry. In practice, however, even a relatively 
small amount of broadening, i.e. F = A/4, is sufficient 
to prevent the correct determination of the OP symme- 
try, as shown in the lower panels of figure [2j unless Z is 
sufficiently large. At low Z, the ZBCP typical of the d- 
wave symmetry, smoothed by F, can be confused with the 
zero-bias maximum of the anisotropic s-wave one. Thus, 
in the presence of a zero-bias maximum in the experi- 
mental curves of afe-plane point contacts in AR regime 
(i.e at low Z) what we can only say is that the OP is zero 
along some directions. Only a&-plane measurements in 
the tunneling regime would allow us to distinguish the 
two situations but, apparently, they are not so easy to ob- 
tain in Fe-based compounds with point-contact or STM 
techniques. 



C. Shape of the Fermi surface and directionality 

In this subsection we will extend the 2D BTK the- 
ory to arbitrary shapes of the FS and will show what 
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is expected when the current is injected either in the ab 
plane or along the c axis. Under a few restrictive hy- 
potheses (still present in this approach), the final result 
will be the full 3D generalization of the BTK model to 
any anisotropic feature both of the FS and of the pair 
potential symmetry. 

The first obvious generalization, essential for the in- 
terpretation of the results in Fe-based superconductors, 
is to consider the multiband nature of superconductiv- 
ity in these compounds and the multiple sheets of the 
FS arising from the crossings of the Fermi energy by the 
different bands. This problem gained for the first time 
importance and popularity ten years ago with the dis- 
covery of MgB2, the first well-known multiband, multi- 
gap superconductor [TBI El- At the level of AR data 
analysis and in order to have a limited number of fit- 
ting parameters the standard approach is to consider a 
two-band, two-gap system where the values of Aj, Zi 
and Fi of every band (i = 1, 2) have to be determined 
by the best fit of theoretical curves to the experimental 
ones. The total normalized conductance is thus written 
as: GtotiE) = WiGi{E) + (1 - Wi)G2{E) where G,{E) 
is the normalized conductance of the i-th band (ID or 
2D depending on the selected model) and wi is the rel- 
ative contribution of band 1 to the total conductance. 
The drawback here is that, apart from a few lucky cases 
[13], the weight factor wi is also a parameter of the fit. 
We will see later that this is not the case in the full 3D 
model where the weighting factors come automatically 
from the complex interaction between the FS shape and 
the direction of current injection. 

A simple extension of the 2D BTK theory to the 3D 
case can be realized under the hypothesis of a spherical 
FS in the S side of the junction (the N side is always 
supposed to have a spherical FS) as shown in [IS]. The 
main point here is to express tn, <Js{E) and A(i5) (if the 
OP is anisotropic) as a function of both the azimuthal 
and inclination angles 9 and 0, respectively. Of course 
the specific expressions depend on the direction of cur- 
rent injection and on the particular symmetry of the pair 
potential. As an example, for a a-axis contact (always 
with no mismatch of Fermi velocity at the interface) it is 
rather easy to show that: 



sph 



(0,0) 



cos^ 9 sin^ ( 



Z2 



and A(i?, 9, <j))± are the expressions of the pair potential 
seen by the ELQs and HLQs, respectively, for the specific 
symmetry considered. By introducing A{E^9^4i)± into 
the above formulas for j± and tpd (equations 6 and 7), 
and then inserting t^^^{9, (/)) and the obtained 7± {E, 9, (p) 
and 'Pd{9, <j>) into equation it is finally possible to inte- 
grate on the two angles arriving at the normalized con- 
ductance: 



sph / 



loll 



The model discussed so far is still implicitly based 
on the assumption of a spherical FS. The problem then 
arises of how to express the c-axis or a-axis conductance 
in a superconductor with a quasi-2D FS. This is the case 
of cuprates or of Sr2Ru04. A possible solution, though 
not really general, is to restrict the integration in </) to a 
narrow angular range around tt/2 as in |15| . In this way, 
the spherical FS of the model is reduced to a almost- 
cylindrical belt that locally approximates the true, cylin- 
drical FS. In the following, we will instead derive and 
present a more general approach to a full 3D case, with 
a FS of arbitrary shape. 

Just for simplicity we will still limit the calculations to 
the contribution of only two bands (thus the band index 
i assumes only the values 1 and 2). As for the direction- 
ality of the contact, we consider here current injections 
along the x axis (j/z-plane interface) and along the z axis 
(xy-plaiae interface), being representative of the typical 
experimental conditions in single crystals of the Fe-based 
compounds, i.e. a6-plane and c-axis contacts, respec- 
tively. From now on the subscript x oi z will denote 
the specific expressions to be used for the two contact 
directions. 

The particular shape of the i-th FS is described by the 
wave vector kp^ i(0, </)), while the unitary vector perpen- 
dicular to the FS at any point of the reciprocal space is 
simply given by np, i(6l, 0) = ^'"'g'^^''^-' x ^^^M^''!'^ whose 
calculation is trivial once the shape of «-th FS has been 
chosen. 

Under the additional (but reasonable) hypothesis that 
the points on the FS are close to the points of maximum 
symmetry of the energy bands (i.e. they are close to 
the top or the bottom of parabolic-like bands) one can 
express the Fermi velocity at any point as a function of 
the component of kp^i perpendicular to the FS at that 
point: 



m 



fikf 



-npa(y, ( 
) • np,i(6 



(8) 



where m* is the effective mass of quasiparticles. Finally, 
the components of the Fermi velocity along a given di- 
rection can be simply obtained by projecting vp_i(0, 0) 
along that direction, i.e. v-px/i{9, (j)) — vp i(0, (/>) • i and 

VFz,i{9,(t)) = vp i(0,(/)) • k where, of course, i and k are 
the unitary vectors of the x and z axis, respectively. We 
have now all the ingredients required to calculate the 
normalized conductance. 

By following the approach originally developed in [HI 
I16j and recently summarized in [S], and by neglecting 
possible interference effects between bands that can lead 
to the formation of bound states at the surface [11] , we 
obtain the normalized conductance for current injection 
along the X axis: 
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where 



[vFx,i{0, 



(10) 



and (ysx,iiE, 9, (p) has the same expression shown m equa- 
tion [5] but now contains the band-associated functions 
l±,i{E,9,(t)), obtained by substituting the OP expres- 
sions A(£^, 6, 4>)±^ i into the 7± expressions (equation 6 



of subsection II B), and the band-associated, direction- 



dependent normal transparency TNi^,i(^,0). 

Equation [lO] directly comes from equation [l] if the dif- 
ference in the effective masses of all the particles and 
quasiparticles involved in the process (in S and N) is dis- 
regarded. The condition of no mismatch of Fermi veloci- 
ties across the interface translates here into wp, 1(6** , 0*) = 
vj^ where 9* and 0* define a point on the FS (usually in 



the k^ky plane) where we suppose the two velocities to 
be equal. By imposing this condition, we eliminate the 
need to know the term h/m* which multiplies all the ele- 
ments at numerator and denominator of equations |9] and 
[To] In this way, of course, we neglect the small deviation 
of the quasiparticles in crossing the interface due to the 
modest mismatch of Fermi velocities in N and S simply 
arisin g fro m the geometry of the FS in S. In equations 



9 and 10 = vn • i = WNCOS^sint/) and Z^^i is the 
barrier height parameter for current injection along the 
X axis in the i-th band. Finally, the limits of integration 
in (j) {(f'min and tt — (pmin) ^rc fixed in such a way that 
the integration is restricted to the first Brillouin zone. 

If all the previous considerations hold, writing the ex- 
pression of the normalized conductance for current injec- 
tion along the z axis is straightforward: 



E,: C '^s.AE, 9, cl>)r^.,d9, 'jyy-^f:$^kl.X9, 0) sin <j>d9dcj^ 



Gmz{E) — 



Er7r/2 ^27r 
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TN2 



(11) 



Here T'^z,i{9,4') is obtained as in equation 10 but by 
projecting onto the k direction. 

It is trivial to show that, in the presence of a single 
band and when not only the FS in N but also the one in 
S is spherical, equations [9| and [TT] reduce to those of the 
3D BTK model already discussed at the beginning of the 
present subsection. 

We can now discuss the results of these calculations in 
a couple of cases of particular interest for Fe-based super- 
conductors. In the typical 1111 or 122 compounds, the 
FS sheets can be approximated by more-or-less warped 
cylinders [T7]. In some special cases, 3D pockets have 
been calculated, that can instead be approximately de- 
scribed as ellipsoids [18]. 

Surfaces that approximate reasonably well these FS 
sheets and that can be easily written in parametric form 
are the oblate or prolate spheroids and the one-sheeted 
hyperboloids of revolution. Here we will use two hyper- 
boloids with different (negative) Gaussian curvatures to 
represent the hole and electron FS sheet of a hypothetic 
1111 or 122 Fe-based compound. These theoretical FS 
sheets are shown as matt surfaces in figure [3];a) and (b). 
As far as the symmetry of OPs is concerned, we will con- 
sider here a isotropic s-wave pair potential on the electron 



FS sheet at M, and a fully anisotropic s-wave symmetry 
with zeros (figure^ (a)) or a dx^-y^-wave symmetry (fig- 
ure [3](b)) for the hole-like FS sheet at F. In figure [s] the 
OPs are represented as semi-transparent gridded surfaces 
with red grids where the pair potential is positive and 
blue ones where it is negative. 

By using these models for the geometry of FS sheets 
and symmetry of OPs and by working out equations [9] 
and[TT]we obtain the theoretical normalized conductance 
for point contacts along the i-axis direction (a6-plane 
contacts) and along the z-axis one (c-axis contacts). This 
conductance is shown as a function of bias voltage at 
T = and for different values of the parameters in pan- 
els (a) and (b) of figure [4] respectively. In all the simu- 
lations we used: Ai = 4 meV (maximum absolute value 
of the anisotropic pair potential), A2 — 9 meV (value 
of the isotropic pair potential), Fi/Ai — r2/A2 = 1/4, 
and ^1 = Z2 = 0.1. In practice the small gap is on 
the hole FS and is anisotropic, with line nodes or zeros 
along four specific directions, while the large one is on 
the electron FS and is isotropic. In addition, Z values 
similar to the experimental ones and reasonably small 
(but not null) broadening parameters are used. The se- 
lected symmetries and parameters do not correspond to 
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FIG. 4: (color online): (a) Normalized AR conductance cal- 
culated by using the full 3D BTK model described in the text 
for the FS shapes and the two OP symmetries shown in Fig. 
3, for two different values of a angle and for current injec- 
tion along the ab-plane; (b) The same as in panel (a) but for 
current injection along c-axis and for two different values of 
Z. The other parameters of the calculations are shown in the 
labels and discussed in detail inside the text. 



FIG. 3: (color online): (a) Reciprocal space representation 
of the shape of the FS sheets (matt surfaces) and of the pair 
potentials (semi-transparent gridded surfaces) used in the cal- 
culations whose results are shown in Fig. 4 (see text). In this 
case the OP on the electron FS sheet at M has isotropic s- 
wave symmetry while the one on the hole FS sheet at F has a 
fully anisotropic s-wave symmetry; (b) The same as in panel 
(a) but for a d^2_y2 symmetry of the OP on the hole FS sheet 

at r. 



any actual Fe-based superconductor. Nevertheless, they 
are somehow inspired by recent theoretical papers on 
hole-doped BaFe2As2 [19] and on BaFe2(Asi_xPx)2 [5^ 
where 3D nodal structures are predicted to appear in 
the warped hole FS. Symmetries involving nodes on the 
electron FS (nodal s-wave) or on the hole FS (rf-wave) 
have been discussed also for 1111 low- Tc compounds 
Figure |4]j a) shows the normalized conductance curves of 
a&-plane contacts for the two selected symmetries of the 
small gap and two different values of the a angle. It is 
clearly seen that, independently of the anisotropic sym- 
metry of the small OP - provided it has zeros in some 
directions - and of the injection angle with respect to 
these directions, the a6-plane curves present a more or 
less pronounced zero-bias peak due to the small gap Ai 
and broader shoulders originated by the isotropic gap 
A2 . This result is somehow expected if one considers the 
curves shown in Fig. 2 for small Z values and F 7^ 0. 
Quite different is the situation of c-axis contacts. As 
shown in Fig. 4 (b) the normalized conductance in this 
case presents no zero-bias peaks, independently of the 
symmetry of the small OP. On the whole, the curves ap- 
pear very similar to what is expected in s-wave symmetry 
but for enhanced Z values (of the order of 0.2-0.4). Also 



the structures due to the small OP have a s-like appear- 
ance in both the symmetries and even for Z = 0, always 
resulting in more or less pronounced shoulders or peaks 
at about 4 meV in the total conductance. This result is 
only partially expected. In fact now the HLQs and the 
ELQs feel the same phase independently of 9 and, thus, 
interference effects due to phase differences cannot occur. 
Nevertheless, the full integration in < < 27r gives rise 
to an average of conductance curves representative of the 
different gap values that are seen at different directions 
in the ab plane. As it happens in the afe-plane contacts of 
panel (a), the contribution of arbitrary small OPs close 
to the nodes or zeros should yield zero-bias maxima in 
the total conductance. But this is not the case. The 
reason is an additional and pure geometric effect. 

When the current is injected along a direction mainly 
parallel to the FS (as in this case, when the current 
is injected along the c axis) the normal transparency 
tnz, K^*, 4') is markedly reduced in amplitude at any angle 
and its overall shape is modified accordingly. It can be 
easily shown that both these effects can be roughly sim- 
ulated by inserting a properly enhanced Z value into the 
TNz{d, (f>) = cos^ (p/ (cos^ (j) + valid for a spherical FS. 
For the FS geometries here discussed this Z-enhancing 
effect occurs also in a6-plane contacts but to a quite mi- 
nor extent since, in this case, most of the FS is almost 
perpendicular to the direction of current injection. In 
other words we can roughly say that, the larger is the 
projection of the FS along the tunneling direction, the 
smaller is the Z-enhancing effect. Of course, as extreme 
limits we have the spherical FS (where the Z enhancing 
is null for every tunneling direction) and the perfectly 
cylindrical FS (along whose axis the enhanced Z tends 
to infinity and, thus, tnz(^, </>) tends to zero). In the lat- 
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ter case, however, the resulting normalized conductance 



is zero as expected from equation 11 and already pointed 
out in [5] . As a consequence of this Z-enhancing effect re- 
lated to the FS-geometry, of the strong Z dependence of 
the zero-bias peaks shown in figure[2]and of the absence of 
phase-difference interference effects, it turns out that c- 
axis contacts on superconductors with warped cylindrical 
FSs cannot show zero-bias conductance peaks or maxima 
even in the presence of d-wave or fully anisotropic s-wave 
symmetry of the OPs. In this case only zero-bias peaks 
or maxima in a6-plane contacts can give information on 
the anisotropic symmetry of one or more of the OPs but, 
as shown in figure |4]^a), without allowing the complete 
determination of this symmetry. In section |III| we will 
provide a couple of examples in which zero-bias maxima 
have been observed experimentally in single crystals and 
polycrystals of Fe-based compounds of different families. 

These results are also interesting from another point 
of view. By reversing the previous reasoning one can 
say that if the OP on a FS sheet has a d-wave or fully 
anisotropic s-wave symmetry and AR experiments in c- 
axis direction do show zero-bias peaks or maxima, it 
means that the Z-enhancing effect is negligible, i.e. this 
FS sheet has a large projected area on the xy plane. In 
materials layered along the ab plane, this occurs when 
the FS sheet has a closed surface, i.e. when it is three- 
dimensional. This situation can be easily simulated by 
substituting one of the two FSs of figure [3] with an oblate 
or prolate spheroid, but a detailed discussion of these 
results is beyond the scope of the present paper. The 
only important message we would like to leave here is 
that in some particular situations, e.g. in presence of 
anisotropic pair potentials, directional point-contact AR 
spectroscopy can provide information not only on the 
symmetry of the OPs but also on the geometry of FS 
sheets where the superconducting gaps open. A less uni- 
vocal and drastic, but still interesting information on the 
FS geometry comes in any case from the observation that 
in the framework of this full 3D BTK model the weight 
factor of every band in the total normalized conductance 
comes directly from the FS topology and the direction- 
ality of the contact and is not a fit parameter as in the 
standard two-band fitting approach. It means that once 
the geometry of FS sheets is fixed (e.g. by first-principle 
band-structure calculations) it is then possible to com- 
pare the experimental PCAR results with the theoretical 
3D BTK predictions getting information on the true FS 
topology as compared to the theoretical one. Just as an 
example, we report here the weight factors of bands 1 and 
2 coming from the results shown in figure|4]and, thus, for 
the geometries shown in figure [Sj They are: = 0.48, 
wf = 0.52, w1 = 0.33 and w§ = 0.67. These numbers 
are not far from the actual values obtained in PCAR ex- 
periments on Fe-based superconductors, particularly as 
far as wf^ is concerned [5DHH]. 

In principle, one should always use this full 3D model in 
fitting the experimental data. However, i) the 3D model 
is rather complicated to handle as a fitting tool; ii) it re- 



quires the knowledge of the FS geometry, which is not al- 
ways available; iii) most of the Fe-based superconductors 
have quasi-2D Fermi surfaces so that the 2D model can 
be safely used as a first approximation. Indeed, the differ- 
ence between the 3D and the 2D models is dramatic only 
along the c axis, and in the case of rf-wave or anisotropic 
s-wave symmetry as shown in figure |4| In the case of 
s-wave symmetry it is possible to show that the results 
obtained with the 2D model are totally compatible with 
those of the full 3D one. For example, we calculated the 
normalized conductances for injection along the ab plane 
and along the c axis by using the full 3D model, using a 
3D FS made of two moderately-warped one-sheeted hy- 
perboloids of revolution. The two s-wave gaps and the 
other parameters used in the calculation are the same 
as in figure [4j except the barrier heights that here are 
Zi = Z2 — 0.2. When the theoretical 3D conductance 
curves are fitted by the 2D BTK model the resulting pa- 
rameters are perfectly coincident with the starting ones 
(including the weights of band 1 and 2) with only the ex- 
pected exception of the parameter Z that results about 
20% larger along ab plane and 2.8-3.6 times larger than 
the actual one in the fit of the c-axis conductance, due to 
the mentioned Z-enhancing effect. As a consequence, in 
s-wave symmetry, the use of the fully 3D model is manda- 
tory only if a precise determination of the Z parameters 
is required. 

In this rather long subsection we have discussed in 
detail the effects of FS geometry and directionality of 
the contact on the results of point-contact AR spec- 
troscopy in multiband multi-gap superconductors like the 
Fe-based ones. In the next subsection we will focus in- 
stead on the effect in PCAR curves of the presence of a 
strong electron-boson interaction and on the role played 
by the shape of the electron-boson spectral function, by 
appropriately solving a three-band s± Eliashberg model 
for the determination of the Ai(_E) functions to be in- 
serted in the above reported expressions for the normal- 
ized conductance. 



D. Electron-boson interaction and shape of the 
spectral function 

It is well known that, in the BCS theory, the OP 
does not depend on energy. This is not the case when 
the Eliashberg theory of superconductivity is adopted. 
By starting from the electron-boson spectral function 
Q!^F(51) it is indeed possible, within the Eliashberg ap- 
proach, to account for the full energy dependence of the 
OP which is now described by a complex function of en- 
ergy, A{E) = ^A{E) + i5A(£'). In particular, as the 
strength of the coupling increases, the imaginary part 
also increases and small deviations from the BCS density 
of states can be observed at the typical boson energies. 
In the moderate/strong-coupling regime these deviations 
can be experimentally detected by tunnelling and PCAR 
spectroscopy. 
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FIG. 5: (color online) (a); Andreev-reflection spectra obtained with different Z values by introducing in the BTK model the 
energy-dependent gaps calculated within a three-band s± Eliashberg model and by using a Lorentzian a^Fiyi) (inset) with 
rJo = 12 meV and half width T = 4 meV. The coupling constant Xtot = 2.11 has been fixed so that Tc = 24.92 K. The gaps 
given by the Eliashberg calculations are Ai = 6.70 meV, A2 — —4.02 meV and A3 — —9.02 meV and correspond to the case of 
Ba(Feo.9Coo.i)2As2 I22j. The other parameters of the BTK model are Zi = 0.3, F^/Ai = 1/4 and Wi = 1/3. The sign-changed 
derivative of the conductance curves are shown in (b). (c): same as in (a) but using a AFSF spectral function (inset) with 
rJo = 12 meV, Xtot = 1.87 and ^max = 4r2o. Calculated Tc is 26.49 K. The gaps given by the Eliashberg calculations are 
Ai = 6.55 meV, A2 = -4.01 meV and A3 = -9.04 meV. The other parameters of the BTK model are Z, = 0.3, T/Ai = 1/4 
and Wi = 1/3. (d): Conductance curves calculated as in (c) but for Z = 0.3 and different values of flmax in the AFSF a'^F{n). 
From bottom to top, Qmas: = 4r2o, ^max = 3^0 and Umax = 2f2o. The three spectral functions are shown in the inset. 



From the theoretical point of view, a precise knovifledge 
of the shape of the spectral function is not really neces- 
sary, at least to a first approximation, to calculate Tc, the 
gaps etc. within the Eliashberg theory, since these quan- 
tities mostly depend on the representative frequency of 
a^F{il) [23;- 25 . However, we will show in the following 
that the shape of the spectral function and the maximum 
boson energy affect the possibility to observe electron- 
boson interaction features in the Andreev-reflection spec- 
tra. In this regard, we will compare different shapes of 
the spectral function, namely a Lorentzian function [26j 
(see inset to figure Isk) and the typical antiferromagnetic 
spin fluctuations (AFSF) spectral function [37] (inset to 
figure [5]^c)). 

The effect of the shape of the spectral function on the 
electron-boson features will be discussed within a three- 
band, s± Eliashberg model which revealed useful in de- 
scribing several Fe-based superconductors [52 [23 E5] . 
Within this model, the electronic structure of Fe-based 
superconductors is described by two hole bands and one 
equivalent electron band in the case of hole doping, and 
by two electronic bands and one equivalent hole band 
in the case of electron doping. The other main assump- 
tions of the model are that: i) the total elcctron-phonon 
coupling constant is small [221 130]; ii) phonons mainly 
provide intraband coupling; iii) spin fluctuations mainly 
provide interband coupling. We will choose, as a rep- 
resentative example, the compound Ba(Feo.9Coo.i)2As2 
[22j . The parameters (A, Z, w etc.) will be indexed 



according to the bands, i.e. 1 for the hole FS, 2 and 3 
for the electron FS sheets. Further details on the model 
and on the parameters of the calculations can be found 
in [22^. 

Figure [sjja) shows several Andreev-reflection conduc- 
tance curves calculated by using a three-band 2D BTK 
model (equation |4]) in which, however, the OPs are not 
constant but energy-dependent. The functions Ai{E) 
were previously calculated within the aforementioned 
three-band Eliashberg model by using coupling constants 
able to give Tc and gaps similar to the experimental ones. 
The a'^F{il) was chosen to have a Lorentzian shape, with 
representative boson frequency fig = 12 meV, similar to 
the spin-resonance energy obtained from neutron scatter- 
ing experiments in Ba(Fei_xCox)2As2 [31]. The curves 
in figure [Sj^a) refer to different values of the Z parameters 
(Zi = Z2 = Z3 ^ Z) from (ideal AR contact with no 
barrier) to 20 (tunnel regime) . Each conductance curve is 
the sum of three contributions (one for each band of the 
model) whose relative weight have been here chosen to be 
equal, wi = W2 = W3 = 1/3. The other parameters of the 
BTK model have been chosen so as to compare to real ex- 
perimental situations, i.e. T = 2 K, Zi = Z2 = .^3 = 0.3 
and Ti/Ai — 1/4. The calculated conductance curves in 
figure |5ja) feature structures at energies higher than the 
gap values, at about Amax + = 21 meV, which be- 
come more pronounced as Z increases. The features can 
be seen much more clearly by looking at the sign-changed 
derivative of the conductance, —dG/dV — —d^I/dV^ 
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shown in figure . A peak is present above 20 meV 
and an additional feature can be also seen around 40 meV 
(vertical dashed lines). When Z decreases, the structures 
slightly reduce in amplitude but their position in energy 
hardly changes, showing that important information con- 
cerning electron-boson interaction can be obtained in the 
Andreev-reflection regime as well. 

Figure [sjc) shows the same cases as in (a) but now cal- 
culated by adopting a normal AFSF spectral function, 
with maximum at flo = 12 meV and maximum boson 
energy {flmax = 4r2o ~ 48 meV). Here the total cou- 
pling constant A is fixed so as to give (approximately) 
the same and the same gaps as in panel (a). The 
additional structures are no longer visible in the con- 
ductance curves at any value of the Z parameter. Only 
tiny deviations from the BCS case can be seen in the 
sign-changed derivative (not shown) but these features 
would be hardly observable in a real experiment. Fig- 
ure [sjd) shows what happens if the value of ^max is 
reduced, though scaling Xtot so as to give the same 
and gaps (see the inset). If ^max = 3r2o small additional 
structures start to become visible and they grow consis- 
tently in amplitude when flmax = 2f2o, but now they 
are at A^^^ -I- ^max and not at A„iax + ^^o- Summariz- 
ing, in this case the electron-boson interaction structures 
are detectable only if the spectral function is abruptly 
cutoff resulting in a shape no more physically plausible 
(see inset to figure [5](d)) and are located at energies that 
have nothing to do with the representative boson energy 
fio- The conclusion is that reasonable parameters for the 
normal AFSF spectral function do not allow reproducing 
features that are actually observed in the experiment (see 
Hence, we suggest that this particular shape 
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of the c?F(J7) is not suitable to catch finer details of the 
PCAR spectra measured in these compounds. This fact, 
however, may not be surprising since the AFSF function 
used here was calculated for a uniform electron gas in 
the normal state [27] which looks too rough an approx- 
imation in this case. Very recently, the AFSF electron- 
boson spectral function in the superconducting state has 
been calculated [35] and is indeed rather similar to the 
Lorentzian a^F(r2) used above, which is thus a reason- 
able and easy-to-handle approximation. 

Figure pfa) shows an Andreev-reflection spectrum 
(solid line) calculated at r = 2 K, Z = 0.3 and Ti/Ai = 
1/4 by using a Lorentzian spectral function, together 
with the partial contributions from the three bands 
present in the model (dash-dot lines). Each of the three 
partial conductance curves shows strong-coupling struc- 
tures at energies higher than the gap. However, while 
in a single-band case one could expect these features to 
appear at A -f JIq (and potentially at higher harmonics), 
in a three-band case with a dominant interband pairing 
the situation is not so straightforward. For instance, the 
contribution of band 2 features a peak at A3 -I- rip ~ 22 
meV and not at A2 -t-^io- This fact, which is most proba- 
bly due to the strong interband character of the coupling, 
has to be investigated further and suggests that quanti- 
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FIG. 6: (color online): Solid line: total conductance curve 
calculated with the same parameters of figure |5|a) but with 
Z=0.3. Dash-dot lines: single band contributions to the total 
conductance. Inset: sign-changed derivatives of the curves 
shown in the main panel, (b) circles: total conductance curve 
from panel (a) compared to the 2-gap 2D BTK fit (solid line) 
with constant BCS values for the gaps. The parameters of 
the fit are reported in the left label. 



tative calculations are recommended in order to reliably 
interpret electron-boson features in Fe-based supercon- 
ductors. 

Figure |6jb) shows what happens when the total con- 
ductance of panel (a), here represented by circles, is fitted 
to the two-band 2D BTK model which adopts constant 
BCS energy gaps (solid line). The agreement between 
the two curves is remarkable in the gap region, despite 
the fact that the BTK model contains only two gaps (a 
three-band model would contain too many parameters 
for the fit to be reliable). The fit allows determining 
very well the small gap (4.0 meV against 4.02 meV, see 
the parameters in the labels) . The large gap obtained by 
the fit is similar to A3 given by the Eliashberg theory, 
with a slight overestimation (9.9 meV against 9.02 meV); 
the effect of the intermediate gap is a hardly visible ex- 
cess of signal with respect to the fitting curve between 
the extremal gaps, which is indeed sometimes observed 
experimentally. At higher energies a considerable devia- 
tion from the fit is observed: The calculated conductance 
shows additional features together with a higher conduc- 
tance up to about 60 meV. Such an excess of conductance 
at intermediate energy has been observed rather often by 
PCAR in Fe-based superconductors (see figure 2d and 2e 
from [33j or figure [Tsjja) below). Furthermore, a higher- 
energy excess conductance that some calculated curves 
show for particular bands can cause an increase of the 
amplitude of the tails even up to 80-100 meV, as shown 
in figure [13)^ a) and in the inset to figure [Tl|a). This effect 
can be a further evidence of the strong-coupling charac- 
ter of these compounds. However, since the height of the 
tails at high energy in the experimental curves depends 
on the normalization procedure, one must be confident 
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that the actual normal state is adopted, i.e. the curve 
at low temperature and H > Hc2- Since the upper crit- 
ical field in these compounds is usually not experimen- 
tally accessible, one could try using the normal state at 
T > Tc, provided that no heating effects are present in 
the point-contact which could cause a downward shift 
of the conductance curves with increasing temperature. 
Therefore, a systematic study of this effect is desirable, 
both from the experimental and theoretical point of view, 
as it could greatly help acquiring more information con- 
cerning the strength of the coupling and the energy of 
the mediating boson in Fe-based superconductors. 

The fact that in the above discussion we have com- 
pared the same data to either a two-band BTK model 
with constant (energy- independent) OPs or a three-band 
BTK model with energy-dependent OPs might look a 
little confusing. Let us just explain this point. 

One one hand, if one starts from an experimental 
PCAR spectrum, the minimal BTK model able to fit 
the data - at least in the small-to-intermediate energy 
range - is the two-band BTK one. In some spectra, a 
discrepancy between the data and the fitting curve is ob- 
served at energies between the small and the large gap, 
which might be due to the presence of a third gap, which 
is however experimentally difficult to discern. Moreover, 
this discrepancy is too small to justify the use of a three- 
band BTK model that would contain 11 free parame- 
ters and, thus, would not be reliable. In this sense, the 
large uncertainty on the gap A2 (especially in the case 
of Sm-1111) might partly arise from the oversimplifica- 
tion of the fitting model. Even if this fact has not pre- 
vented obtaining interesting information by PCAR mea- 
surements, this consideration has to be kept well in mind 
when studying these complex multiband systems. 

On the other hand, the minimal Eliashberg s± model 
that allows obtaining the experimental gaps and Tc is a 
three-band, three-gap one as explained in |26j . If one 
wants to take into account the effects of the energy de- 
pendence of the OPs and of the strong electron-boson 
coupling on the shape of the PCAR spectrum, one is 
thus forced to generate the Andreev-reflection curve by 
using a three-band BTK model with the three OPs com- 
ing from Eliashberg calculations. Note that, in the latter 
case, the OPs are not fitting parameters and only F^, Zi 
and the weights can be adjusted to "fit" the experimental 
spectrum. 



III. SYMMETRY OF THE ORDER 
PARAMETERS IN FE-BASED 
SUPERCONDUCTORS 

A. "1111" compounds 

In most of the 1111 compounds, the model for spin- 
fluctuation mediated superconductivity predicts a s± 
gap symmetry with a sign change between hole-like and 
electron-like Fermi surface sheets. This is indeed com- 



patible with what has been found by point-contact mea- 
surements, despite initial contradicting claims. An order 
parameter with nodes on the electron FS has been pre- 
dicted in 1111 compounds with low Tc [5J which have not 
been studied yet by PCAR. 

The first PCAR measurements in 1111 materials like 
LaFeAs(0,F) and SmFeAs(0,F) showed zero-bias peaks 
that were initially interpreted as being due to a nodal 
gap symmetry, namely the d-wave one. However, the ap- 
pearance of these zero-bias features and their height were 
later discovered to depend on the pressure applied to the 
sample by the tip (in the needle-anvil configuration) [33] ; 
in some samples, they were frequently found in some re- 
gions of the sample surface and never in others; moreover, 
their temperature dependence was odd and completely 
inconsistent with that expected for the zero-bias peak 
associated with Andreev bound states in nodal supercon- 
ductors [35| I36j . At the same time and most importantly, 
increasing pieces of evidence were collected that, in suit- 
able conditions (better sample quality, or zero pressure 
applied to the sample as in the "soft" technique described 
in |9j) the zero-bias peak did not show up ^2ni [HI [55] . 
In the end, the majority of the scientific community got 
convinced that the zero-bias peaks early observed in 1111 
compounds were not a signature of d-wave symmetry of 
the order parameter. We will thus not discuss these re- 
sults any longer; a short review on this subject can be 
found in [5]. 

The absence of zero-energy bound states (particularly 
in a&-plane contacts) certainly excludes those symmetries 
where nodal lines lie on the Fermi surface. However, the 
predicted symmetry of the OP in pnictides does involve 
a change of sign in the OP, but this should occur be- 
tween hole-like and electron- like Fermi surface sheets [3] . 
According to some models, this sign change should give 
rise to observable effects in the PCAR spectra through 
interband interference effects. The article by Greene et 
al. in this same issue is specifically devoted to exploring 
these models and comparing their predictions to experi- 
ments. Here we will not deal with this aspect and show 
instead how the most reliable PCAR spectra in 1111 can 
be successfully fitted by simpler models that disregard 
interference effects. Indeed, the predicted features as- 
sociated with interband interference (i.e. surface bound 
states, both at zero and nonzero energies [TT]) may be of- 
ten made indiscernible by various sources of broadening 
or overwhelmed by stronger effects. 

It is however worth mentioning that, in contacts with 
vanishing barrier at the interface (let's say, Z = in the 
BTK language), a strong reduction in the Andreev sig- 
nal is predicted [TT] as a result of destructive interference. 
In Fe-based superconductors a smaller amplitude of the 
Andreev signal than in more conventional systems (even 
multiband, like MgB2) is indeed often (but not always) 
observed experimentally. However, a similar reduction in 
the amplitude of the superconducting signal is also ex- 
pected if the normal bank of the contact is a diffusive 
metal [371 [3H] , if a large inelastic scattering occurs at the 
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interface [35], if a fraction of the probe current is injected 
in normal regions (as it can happen if the dopant distri- 
bution is inhomogeneous) , or even if some Fermi surface 
sheets are gapless. For the sake of completeness, it must 
be said that a small Andreev signal is also predicted if 
the gap is anisotropic (see figure [2| . At this stage of the 
research it is impossible to say which of these effects is 
really taking place. This is particularly true for 1111, 
since the growth of large single crystals is still very diffi- 
cult and most of the experiments have been carried out 
on polycrystalline samples which necessarily do not al- 
low ruling out extrinsic effects due to grain boundaries, 
surface degradation, inhomogeneous doping and so on. 



1. LaFeAs(Oi-^F^). 

LaFeAs(Oi_a:Fa:) is chronologically the first Fe-based 
superconductor studied by point-contact spectroscopy, 
and also the first on which we performed "soft" PCAR 
measurement (where the contact is made by placing a 
small drop of Ag paste on the sample surface instead of 
pressing a sharp metallic tip against it) to avoid pressure 
effects that could possibly give rise to spurious zero-bias 
peaks, as shown in some early measurements reported 
in literature. 

We studied polycrystalline samples with nominal F 
content x = 0.1 and thus at optimal doping. However, 
the local F concentration in the crystallites was found 
to show deviations 6x = 0.02 around this value. Indeed, 
point contacts on different regions gave critical tempera- 
tures (defined as the temperature where the Andreev- 
reflection features disappear and the normal-state spec- 
trum is recovered) ranging from 27.3 K to 31.0 K. These 
temperatures lie in the upper part of the resistive tran- 
sition (see figure [rj^e)). 

In most of the contacts we observed clear conductance 
maxima, symmetric about zero bias and related to a 
small gap Ai, together with additional structures that 
we interpreted as being due to a larger order parameter 
A2 (although we acknowledged the possibility that they 
could be features of the normal state, which is also chang- 
ing with temperature [5DJI3D]). In some contacts, these 
additional features were extremely clear and sharp (see 
figure 2 in [20]); in others (see figure [tJ^ a)) they showed 
up as broad shoulders progressively decreasing in height 
and width on increasing temperature. Other peculiar 
features of the spectra are the right-left asymmetry, and 
the residual upward curvature in the normal-state con- 
ductance that looks like a "pseudogap" feature and is 
progressively reduced on increasing the temperature. Cu- 
riously enough, the temperature at which it disappears, 
leaving a flat but asymmetric conductance, is close to 
the temperature of magnetic ordering in the parent com- 
pound. 

The origin of this behaviour is still not clear and will be 
possibly clarified only when single crystals are available 



(see also III A 2 for further discussion) . 



The complete absence of zero-bias peaks or maxima 
such as those shown in figure |4] in the conductance, to- 
gether with the random orientation of the grains, made 
us conclude that none of the order parameters shows line 
nodes or zeros on the Fermi surface. The further conclu- 
sion in favor of a isotropic gap symmetry came instead 
from the fit of the conductance curves, which was indeed 
always possible by using two isotropic order parameters. 

The fit requires of course a normalization of the con- 
ductance curves, which is particularly crucial in these 
materials (as already mentioned in section II D I because 
the gap or boson structures can extend to very high en- 
ergy. In some spectra, these structures become clearly 
visible only if a proper normalization is chosen. In prin- 
ciple, the normalization should be obtained by dividing 
the superconducting conductance curve measured at a 
given T < by the normal-state one at the same tem- 
perature. The latter is experimentally inaccessible due to 
the large critical fields of these materials [32]. Some au- 
thors then choose to divide the raw conductance curve at 
a given temperature by a straight line ^33; that connects 
its high-energy "tails". This procedure generally allows 
a good estimation of the small gap alone. As a mat- 
ter of fact, the pseudogaplike curvature of the normal- 
state conductance curves observed at T > is presum- 
ably present also below (and likely to become more 
and more pronounced when the temperature is lowered) . 
This curvature is completely disregarded by the afore- 
mentioned normalization procedure. To partially correct 
for this, a possibility is to divide all the raw conductance 
curves in the superconducting state by the normal-state 
curve measured just above ^ or (as we also did in La- 
1111) to try to mimic the temperature dependence of the 
(presumed) normal state [20]. We showed that choosing 
one or another of these two procedures changes very lit- 
tle the value of Ai extracted from the fit (less than 2%), 
while the value of A2 can change by about 10%, though 
preserving its trend as a function of temperature, mag- 
netic field and critical temperature. 

Figure [7][c) shows the low-temperature experimental 
curve of panel (a) after normalization (i.e. division by the 
normal state at T^), together with its best-fitting curve 
calculated within the two-gap 2D BTK model discussed 
in section [II A| Let us just recall here that this model is 
really suited to describe a 3D superconductor only if the 
gap is isotropic and the Fermi surface is spherical, which 
is clearly not the case in these materials [ITj [43] . More 
sophisticated models accounting for the real shape of the 



Fermi surface (see section II C I should thus be used, but 



calculations may become too heavy for a fitting proce- 
dure to converge. In practice, using an approximated 
2D model often simply results in an overestimation of 
the barrier parameter Z, as already mentioned in section 
III CI The values of the OPs extracted from the fit of the 
low-temperature conductance curve are Ai = 3.8 ±0.1 
meV and A2 = 10.0 ± 0.2 meV. The fit of another set of 
curves measured in a point contact with = 28.6 K, 
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FIG. 7: (color online): Point-contact spectra in LaFeAsOi-^Fi polycrystals [20] and films [41 j . (a, b): raw conductance curves 
of two contacts with similar critical temperature obtained in polycrystal and film, respectively. In (b), the normal-state 
conductance was not measured and is simulated by the bottom dashed line, (c, d): the low-temperature conductance curves 
from (a) and (b) after normalization (i.e. division by the normal-state curve at T'^\ and the relevant two-gap 2D fit (lines), 
(e): comparison of the resistivity of the polycrystal (left axis) to the resistance of the film (right axis), (f): Temperature 
dependence of the gaps as obtained from the 2-gap 2D fit of the curves in (a) (■) and from the fit of the PCAR spectra of a 
different contact, with — 28.6 K (□), on the same sample. The results of the 1-gap BTK fit of the curves in (b) performed 
by Naidyuk et al. is also shown (o) as well as the 2-gap fit of their low-T curve shown in (d) (T, A). The lines are BCS-like 
curves that fit the low-temperature trend of the gaps and, once extrapolated, would give an "apparent" critical temperature 
T* < T^. 



but featuring clearer structures related to the second gap, 
gives instead Ai = 2.75 ±0.1 meV and A2 = 7.9 ± 0.2 
meV. How can these result be reconciled with the s± 
picture? As already mentioned above (section II D ) we 



have shown that the s± picture is compatible with 
the opening of gaps of different amplitude on the two 
hole-like FS sheets and on the electron-like sheet of hole- 
doped 1111 compounds. Eliashberg calculations assum- 
ing a mainly-interband coupling due to spin fluctations 
give, for the case of LaFeAs(0,F) analyzed here, three 
isotropic order parameters of different signs: a small pos- 
itive OP Aa — 2.82 meV (presumably on the outer hole 
FS cylinder around F), a large positive OP Af, — 8.01 
meV (on the inner hole FS cylinder), and a large neg- 
ative OP Ac = —7.71 mcV (on the electron FS sheets 
around M) [35] . In the absence of interband interference 
effects, Andreev reflection is only sensitive to the ampli- 
tude (i.e. the absolute value) of the OPs (which is the 
actual energy gap). This, together with the similarity be- 
tween the absolute values of the two large OPs, reconciles 
the findings of PCAR with the theoretical expectations 
in the s± model. 

The gaps (or, better, the relevant 2/S./kBT^ ratios) 
extracted from the fit of the conductance curves in fig- 
ure [7|[a) are shown as solid squares in panel (f ) . Their 
temperature dependence is anomalous; both Ai and A2 
follow a BCS-like trend up to T* ~ 18 K, but above this 
temperature A2 is no longer detectable and Ai shows a 
"tail" ending up at the real critical temperature of the 



contact, = 27.3 K. The same result was found in 
other contacts with different critical temperature (20j . 
The data shown by open squares in[7][f) refer to a con- 
tact with — 28.6 K on the same sample; also in this 
case the large gap tends to zero at T* ~ 21 K, and above 
this temperature the small gap shows a "tail" up to the 
real T^. 

Very similar results have been recently obtained by 
Naidyuk et al. [H] who performed PCAR spectroscopy 
with a Cu tip in LaFeAsOi-j^F^: films with Tc = 28 K 
(defined as the temperature where the resistance drops to 
90% of its normal-state value, see fLgure[7][e)) and a tran- 
sition width 5Tc ss 6 K. The conductance curves dl/dV 
obtained from their data are shown in figure [7]jb) , and 
show two clear maxima related to a small gap Ai. The 
authors normalized the curves to a parabolic curve that 
fits the data for \V\ > 10 mV. This necessarily makes all 
the features at \V\ > 10 become small oscillations around 
1 in the normalized curve. As a matter of fact, the curves 
normalized in this way can be fitted by a single-band, s- 
wave BTK model including a broadening term and the fit 
thus gives only a small gap, shown in panel (f) as circles. 
The gap follows a BCS-like curve that, once extrapolated, 
would give again a critical temperature T* ~ 18 K. Struc- 
tures in the raw conductance curves, although impossi- 
ble to fit with the same model, persisted well above that 
temperature and at least up to 25 K, as shown in panel 
b. As a matter of fact, the temperature dependence of 
the raw conductance curves is not dissimilar from that 



15 



observed by us in polycrystals (panel a). Let us show 
what happens instead if one normaUzes the same curves 
by the normal state at . Although the corresponding 
curve is not reported in [41| , based on the similarity with 
panel (a) one can guess that it should look like the dashed 
curve in panel (b). If this is used to normalize the low- 
temperature spectrum, one obtains the curve depicted 
in panel (d), which shows a rather high Andreev signal 
and shows clear shoulders in addition to the structures 
related to Ai. The two-band 2D model allows fitting 
this curve very well, giving Ai ~ 2.9 mcV and A2 ~ 9.5 
meV (triangles in panel f). The agreement between the 
data sets measured by us in polycrystals and by Naidyuk 
et al. in films is incredibly good if the raw conductance 
curves are treated in the same way. This proves that: i) 
the different sets of data are really compatible with each 
other; ii) even if in the raw data by Naidyuk et al. the 
large-gap features are not as clear as in ours, a suitable 
normalization makes them show up clearly enough to be 
fitted by the 2D model; iii) the gap values are rather 
robust against the sample form and synthesis technique, 
provided that the critical temperature is similar, as in 
this case. The two data sets, however, also pose some 
questions on whether the observed anomalous tendency 
of the gaps to close around T* is an artifact or is intrin- 
sic to this material. The simplest possibility, of course, is 
that the behavior of the gap is the result of the shorten- 
ing of the mean free path on increasing temperature, so 
that the junction ceases to be ballistic at a voltage that 
decreases with temperature 



2. SmFeAs(Oi-^F^). 

Also in Sm-1111 the PCAR results reported in litera- 
ture apparently give contradicting results, especially as 
far as the observability of multiple gaps is concerned. 
Here we will show that, however, raw data measured in 
different kinds of samples by different groups can give sur- 
prisingly consistent results if they are treated in the same 
way. Figure[8](a) shows the raw conductance curves mea- 
sured by using the "soft" PCAR technique in a optimally 
doped SmFeAsOo.8Fo.2 polycrystal with Tc — 52 K. We 
can notice that: i) the low temperature curve shows max- 
ima at about ±5 mV, plus additional structures (here in 
the form of shoulders) at about ±15 mV; ii) the normal- 
state conductance (thick line, corresponding to T = 52.1 
K) is not flat but shows a residual, asymmetric "hump" 
around zero bias, which is progressively smoothed on in- 
creasing temperature. In all the soft point contacts we 
made in these polycrystals, this hump disappears - leav- 
ing a flat, though still asymmetric conductance - at about 
140 -150 K (which is the temperature of magnetic order- 
ing in the parent compound) [HI |33] . The origin of this 
"crossover" is unclear. One could doubt that it might 
be related to some residual magnetic order, but no trace 
of a secondary magnetic phase with T^v ~ 140 K has 
been observed in our samples at optimal doping (see e.g. 



|45)). A structural transition has been claimed to occur 
between 175 K and 125 K along the whole phase diagram 
of Sm-1111 [33], but we don't have evidence of this effect 
in the samples studied here. A small anomalous kink in 
the resistivity at similar temperatures, common to fill 
compounds (see e.g. |47]), has been instead observed, but 
its origin is not well understood yet. Transport measure- 
ments 4Sj and ARPES [49^ in detwinned BaFe2As2 crys- 
tals showed that in-plane electronic anisotropy occurs at 
temperatures above the structural transition, thus indi- 
cating the presence of fluctuation effects [15] which have 
been interpreted as spin anisotropy in the paramagnetic 
phase [SD]. It might also be that similar effects are play- 
ing a role in the 1111 systems and are possibly connected 
to the kink in the resistivity and to the shape of the point- 
contact conductance. However, detailed studies have to 
be performed as soon as bigger crystals are available to 
clarify these points. Generally speaking, we can say that 
the proximity of superconductivity to the AF SDW state 
of the parent compound suggests the persistence, in the 
doped compound, of spin fluctuations which could possi- 
bly give rise to the observed behavior of the normal state 
at high temperatures. However, since there are no clear 
predictions concerning the expected shape of the PCAR 
conductance spectra and of their evolution with tempera- 
ture in the situations reported above, further experimen- 
tal and theoretical investigation should be performed in 
order to clarify this issue. Then, it could also be possible 
to better understand whether the different shape of the 
normal-state conductance in La-1111 and Sm-1111 poly- 
crystals is due to different causes or it is only a different 
appearance of the same phenomenon. 

Anyway, the anomalous shape of the normal-state con- 
ductance makes again the normalization be a problem. 
Dividing by the curve at is here a rougher approxi- 
mation than in La-1111 because the critical temperature 
is higher and one could reasonably suspect that the nor- 
mal state at can differ very much from that at low 
temperature. The result of such a normalization is shown 
in panel (c) for the curve at the lowest temperature. The 
normalized curve still presents a strong right-left asym- 
metry, that some authors remove by a symmetrization 
procedure. We chose instead to separately fit the left 
and right side of the curves with the two-band 2D BTK 
model, as shown by the lines superimposed to the sym- 
bols in panel (c). The values of the small gap Ai ex- 
tracted from the two fits (5.7 meV) coincide, but the 
values of A2 are rather different: 17.0 and 19.0 meV, 
respectively. 

The fit of the symmetrized conductance obtained from 
the one reported in figure [sj^c) gives Ai = 5.7 meV and 
A2 = 18 meV. While the small gap remains the same the 
large gap is actually an average of those obtained from 
the separate fit of the right and left side of the unsym- 
metrized one. The symmetrization is thus an acceptable 
practice when the asymmetry of the conductance curve is 
not too large. This is what happened in case of Co-doped 
Ba-122 (see [22 J where we symmetrized the conductance 
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FIG. 8; (color online): PCAR spectra obtained in SmFeAs0i_2,Fi polycrystals [2T]|42(a) and films [4T](b). The bottom curve 
in (b) is an estimation of the normal state conductance. (c,d): the low-temperature conductance curves from (a) and (b) 
normalized to the normal-state conductance at and the relevant two-gap, 2D BTK fit. In (c) the left and right side have 
been fitted separately; in (d) the curve has been symmetrized, (e): Normalized gap values (2A/fc_gr^) as a function of T /T^ 
as obtained in different samples: optimally doped polycrystals (A,V, from ^21,); underdoped polcrystals (o from [3T], -A' from 
[33]); single crystals (DjO from |5l]); films (• from [41]). The green stars indicate the results of the fit shown in panel (d). 
Lines are BCS-like curves. The dashed line fits the trend of the gap taken from [?T], which closes at the "apparent" critical 
temperature T* < . Colored regions approximately indicate the region where the values of Ai and A2 are spread. 



spectra. Since m our results on Sni-1111 the curves show 
a rather large asymmetry we consider more reliable, in 
that specific situation, a separate fit. 

The temperature dependence of the OPs (normalized 
as in figure [7|[f)) obtained in this way is shown in panel 
(e) by up (left side) and down (right side) triangles. Un- 
like in La-1111, the gaps follow a BCS-like trend up to 
T^, that coincides with the bulk EJ- The same panel 
also shows the gaps obtained by fitting the negative-bias 
side of conductance curves obtained in a different poly- 
crystal, with a different doping content (x = 0.09) and 
thus a smaller = 39 K (open circles). Despite these 
differences, once the temperature and energy scales are 
normalized as in panel (e), the values of the gaps agree 
very well with those obtained for x = 0.20 [Uj. A 
few measurements of soft PCAR spectroscopy were also 
performed in very small Sm-1111 single crystals (about 
200 X 200 ^m^) with the current (mainly) injected along 
the ab planes [ST]. As usual, we fitted the negative-bias 
and the positive-bias side separately to account for the 
observed strong asymmetry of the spectra. At low tem- 
perature, the gaps turned out to be Ai = 6.45 ± 0.25 
meV and A2 = 16.6 ±1.6 meV, while = 50.1 K. The 
normalized values of the gaps (actually, the values of the 
gap ratio 2A/kBT^) are shown as a function of T/T^ in 
figure [sf^e) as squares and diamonds. The agreement of 
the data with those obtained in polycrystals is remark- 
able. In general, when all the data sets obtained in var- 
ious polycrystals and crystals are considered together, 
the values of the 2Ai/kBT^ agree very well with one 
another, varying only between 2.5 and 3 (i.e. within the 



thinner colored region in panel (e)). As for 2A2/kBT^, it 
is more scattered and varies at low temperature between 
7 and 9 (i.e. within the wider colored region). This larger 
spread arises both from the asymmetry of the curves and 
from the greater uncertainty in the determination of the 
large gap, also related to the uncertainty introduced by 
the normalization or to the possible presence of two large 
gaps of similar amplitude. 

PCAR measurements in the needle-anvil configuration 
were performed by Chen et al. in polycrystals with a re- 
duced Tc = 43 K with respect to the optimal-doping one 
[551 [S2]. Most of the low-temperature curves shown in 
[33] resemble rather closely those measured by us; they 
present clear peaks related to a small gap plus addi- 
tional shoulders and structures. Unlike in our measure- 
ments, the normal-state conductance above presents 
a pseudogaplike upward curvature similar to that re- 
ported above for La-1111. Prior to fitting their curves, 
the authors use a normalization procedure that is equiva- 
lent to supposing that the normal state is a straight line 
fitting the tails of the unnormalized curve. The prob- 
lem is that, the higher is the temperature, the smaller 
is the voltage range where this approximation holds (see 
Supplementary Information in [33j ) . Moreover, this obvi- 
ously implies disregarding possible structures at energies 
much higher than that of the small gap. Traces of the 
(presumable) large gap remain in the curves, but become 
so small to be confused with minor features, possibly re- 
lated to the electron-boson coupling (sec II D and IV I. 
Instead, if the curves are divided by the normal-state 
conductance, the features connected to the large gap are 
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enhanced with respect to other features at even higher 
energy. The fit of the conductance curves performed by 
Chen et al. gives indeed a single gap, which follows an 
almost perfect BCS temperature dependence shown in 
figure |8][e) as solid stars (note that, because of the differ- 
ence in critical temperature between different data sets, 
the axes are normalized) . In |52| the same authors notice 
that the energy of some extra features (which at T = 4.5 
K occur at about 20, 40 and 55 meV) evolve in a BCS- 
like way on increasing temperature so that, if properly 
normalized to its low-T values, all the curves are super- 
imposed to the gap. On the basis of what discussed in 
section IIP it is tempting to interpret these features as 



is indeed suggestive of strong electron-boson structures 
in the conductance that could only taken into account 
by inserting the energy-dependent order parameters into 
the 2D BTK model (see section 



II D) 



the hallmarks of the strong electron-boson coupling. 

More recently, PCAR measurements in the needle- 
anvil configuration have been performed in Sm-1111 films 
by Naidyuk et al. [41]. Figure |8]^b) shows the tempera- 
ture dependence of one of their curves, which show very 
clear structures related to the small gap and a pseudo- 
gaplike normal state. The peculiar evolution of the tails 
of the curves as a function of temperature seems to in- 
dicate a large excess conductance at high energy, possi- 
bly due to electron-boson structures, as discussed in sec- 
tion |IID Features connected to the superconductivity 
are clearly present up to ~ 30 K. This is the tem- 
perature where the resistance has dropped to about 80% 
of the normal-state value, although the transition starts 
at about 34 K. As in La-1111, the authors' normaliza- 
tion makes the Andreev signal decay to 1 already at ±12 
meV, thus flattening all the structures at higher energy. 
This, again, allows a fit of the small-gap features alone, 
and only to about 18 K; above this temperature the nor- 
malization produces structures that cannot be fitted by 
a BTK model. The gap values follow a BCS-like trend 
that, if extrapolated, would give a critical temperature 
T* ~ 21.6 K. This, together with the low-temperature 
value A = 3.35 meV, gives an almost perfect BCS ratio 
2A/A;bT* = 3.6. However, the "true" critical tempera- 
ture of the contact is ~ 30 K, that gives instead a 
gap ratio of only 2.6, which is exactly the same we ob- 
tained for the small gap (see figure [Sjje)). Based on this 
similarity, we tried to normalize the curves in |41j follow- 
ing our method, and then divided all the curves by the 
presumable normal state at (dashed curve in panel 
b). The result of this normalization (and of a subse- 
quent symmetrization) is shown by symbols in panel (d) . 
The curve has then been fitted to a single-band 2D BTK 
model (dashed line) and to a 2-band 2D model (solid red 
line). It is clear that the two-gap model allows a better 
fit of the small-gap feature's width but also catches the 
position of the abrupt change of curvature, which turns 
out to be related to the second large gap. The values of 
the gaps given by this tentative fit are Ai — 3.5 meV 
and A2 = 11 meV. Their values, properly referred to the 
relevant T^, are shown in panel (e) as green stars and 
again turn out to be in very good agreement with ours. 
The failure of the 2D BTK model to fit the tails of the 
curve in (d) is not too worrying; the shape of these tails 



B. "122" compounds 

In 122 compounds, the model for spin-fluctuation- 
mediated superconductivity predicts a nodeless OP 
which, according to some authors, can however evolve in 
a peculiar nodal symmetry, with three-dimensional nodes 
on one hole-like Fermi surface, when the latter acquires a 
more three-dimensional character [5 . In a recent paper, a 
relationship is claimed between 2Ah/kBTc (being Ah the 
gap on the hole-FS) and the occurrence of nodes in the 
electron FS [53]. According to these authors, both the 
hole-doped system (Ba,K)Fe2As2 and the electron-doped 
system Ba(Fe,Co)2As2 present nodeless gaps around op- 
timal doping, while a nodal OP is very likely in the 
isovalent-doped compound BaFe2(As,P)2 [51 153). 

122 compounds are more suited to PCAR measure- 
ments than 1111 compounds, since they can be grown 
in the form of large single crystals that can be cleaved 
more or less easily. This allows directional measurements 
with the probe current injected either in the basal plane 
or along the c axis. Moreover, in some of these com- 
pounds, ARPES measurements have provided valuable 
information on the amplitude and symmetry of the gap 
on the various sheets of the FS. This information gener- 
ally confirms or supports the finding of PCAR in these 
materials. 



1. (Ba,K)Fe2As2 

ARPES experiments in the hole-doped system 
(Ba,K)Fe2As2 at a; = 0.4 (T^ = 37 K) [Ml [551 showed 
the presence of four isotropic gaps: three of similar am- 
plitude (between 11 and 13 meV) on the inner hole- like 
FS cylinder and on the two electron-FS sheets around 
M, plus a small one (~ 5.8 meV) on the outer hole FS 
sheet. Within the uncertainty of the measurement, which 
is rather large, the gap amplitudes were shown to be 
compatible with a symmetry A(fc) = Aq cos(A:a;) cos(fcj,) 
which is indeed the functional form of the s± symmetry 
in the reciprocal space. The gaps seemed not to follow a 
BCS-like curve as a function of temperature. 

The first directional PCAR measurements in the same 
system were carried out by Samuely et al. |561 15 7j who 
used single crystals of Bao.55Ko.45Fe2As2 with supercon- 
ducting transition between T°" = 30 K and T^''''° = 27 
K. By injecting the current along the ab planes, they 
obtained spectra with very clear two-gap structures and 
no trace of zero-bias peaks or maxima, thus ruling out 
the possibility of any symmetry of the OP involving a 
change of sign or zeros on a single FS sheet. The conduc- 
tance curves were systematically normalized to the con- 
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ductance curve at T^, which sometimes shows a hump- 
like structure (similar to that observed by us in Sm- 
1111 polycrystals) and sometimes a pseudogaplikc, V- 
shaped feature. For current injection along the c axis, 
a similar pseudogaplikc feature and no Andreev signal 
were observed. The zero-bias conductance minimum was 
shown to be progressively filled on increasing tempera- 
ture, up to the temperature (about 85 K in these sam- 
ples) where weak signs of the onset of a magnetic or- 
der are found in resistivity and specific heat |56) . The 
most reasonable interpretation of these findings is that 
magnetic order and superconductivity coexist and are 
probed sometimes together, sometimes exclusively, by 
PCAR [56]. This possibility is confirmed by direct ev- 
idences of such a phase separation on a scale of about 
50 nm [58J, and also by successive PCAR measurements 
in Bao.6oKo.4oFe2As2 crystals with Tc ~ 37 K carried 
out by Lu et al. j59]. These authors indeed observed 
a typical V-shaped conductance valley with a universal 
functional form G{V) = G(0) + c|y|" being n ~ 2/3 
in a fraction of their c-axis contacts on this compound, 
but also in the nonsuperconducting parent compound 
BaFe2As2, as weU as in (Sro.6Nao.4)Fe2As2 (Tc = 36 K) 
and Sr(Feo.9Coo.i)2As2 (T^ ~ 22 K). 

Let us go back to the a6-plane measurements in 
(Ba,K)Fe2As2 by Samuely et al. After the normaliza- 
tion, the conductance curves were fitted to a two-band 
2D BTK model with two isotropic gaps. In different con- 
tacts, the small gap Ai was found to vary between 2 and 
5 meV, and the large gap A2 between 9 and 10 meV. 
The temperature dependence of the gaps was found to 
agree very well with BCS-like curves. The beautiful spec- 
trum shown in figure |9|a) gives 2Ai/A:bT^ ~ 2.7 and 
2IS.2/kBT^ ~ 9, in very good agreement with the find- 
ings in La-1111 and Sm-1111 shown in figures [7][f) and 
[sjje). In some cases, small-amplitude curves with a sin- 
gle broad maximum were obtained, that can be fitted by 
the same model with very large values of the broaden- 
ing parameters Fi and F2. Let us note that these lat- 
ter spectra could be mistaken for evidences of nodal or 
anisotropic OPs if there weren't clearer spectra ruling 
out this possibility. This indicates how the cleanliness of 
the surface is essential for the observation of clear gap 
structures (and thus a correct interpretation of the data) 
as also pointed out by Lu et al. [SS], who performed c- 
axis PCAR measurements in Bao.6oKo.4oFe2As2 crystals 
either cleaved or uncleaved. In uncleaved crystals, they 
observed spectra with a single, ill-defined maximum at 
zero bias and shallow structures at about ±15 meV. In- 
stead, in freshly cleaved samples two clear conductance 
maxima were observed in the spectra. Once normal- 
ized to the normal-state conductance (showing a hump- 
like structure) these spectra were fitted to a single-gap 
BTK model giving Ai = 3.0 — 4.0 meV, corresponding 
again to 2Ai/fcsT'^ — 2.6 and thus in good agreement 
with Samuely's measurements along the ah plane. In ad- 
dition to the clear small-gap features, the raw spectra 
also presented smaller structures that, in the normal- 



ized curve, show up as shoulders of very small ampli- 
tude. These structures were interpreted by Lu et al. as 
being "dips" due to the non-perfectly ballistic conduc- 
tion regime through the contact, in turn due to the very 
small electron mean free path in these compounds |59j . 
They thus argued that c axis PCAR measurements do 
not allow to detect the large gap because it pertains to 
cylindrical FS sheets (the inner-holelike at F and the two 
electron-like at M), while the small gap opens up on a 
more 3D Fermi surface sheet [60] . 

In figure [9] we compare two representative spectra 
taken along the ah plane [SH] H7| and along the c axis 
[59j , by plotting them as a function of the normalized en- 
ergy eV/ksTc (this is necessary because of the different 
Tc's of the samples). The similarity between the spec- 
tra is striking; not only the small-gap features, but also 
the additional structures now lie in the same position. 
This suggests that the shallow structures at ±15 meV ob- 
served in the c-axis spectrum might not be an artifact (as 
a matter of fact, they always occur at the same voltage 
irrespective of the junction resistance |59j) but might in- 
stead be related to the large gap. If this is the case, their 
small amplitude may be explained by the much smaller 
"weight" of the relevant bands than in the ah plane (see 
section II C I . The normalized curves are reported (again 
as a function of eV/ksTc) in figure[9j^b) as open symbols. 
The upper curve is compared to the 2-band fit by Szabo 
et al. [57 , and the relevant gap amplitude and weight 
are reported. The lower curve is instead compared with 
a two-band 2D BTK fit that gives the same small gap 
obtained by Lu et al. (Ai =4.1 meV) but also a large 
gap A2 ~ 16 meV. Note that the weight of the bands 
featuring the large gap is now only 0.1. Although this 
value of A2 seems huge, the ratio 2A2/A:bT^ is the same 
as for the a6-plane curve. If our interpretation is cor- 
rect, the available PCAR experiments in (Ba.K)Fe2As2 
thus concur to a picture of multiple nodeless gaps, in 
agreement with the findings of ARPES. They also seem 
to indicate that these gaps are isotropic (at least, this is 
what one would infer from the directional invariance of 
the gap ratios 2Ai_2/fcsr^^). 



2. Ba(Fei-^Co^)2As2 

The first results of PCAR measurements in 
Ba(Feo.93Coo.o7)2As2 single crystals with Tc — 23 
K are reported by Samuely et al. 56J. The c-axis 
spectra show clear Andreev signals, and in one case 
a double-peak structure, but no evidence of multiple 
gaps. They were thus fitted to a single-gap, s-wave 
BTK model giving A = 5 — 6 meV which corresponds 
to 2A/kBT^ = 5.3 - 6.3 since = 22 K. A large 
broadening parameter (F > 0.5A) was necessary to fit 
the curves. An example of these curves is shown in figure 
^c) (curves 1 and 2). In [59], also Lu et al. reported 
on c-axis PCAR measurements on Ba(Feo.9Coo.i)2As2 
crystals. They observed V-shaped conductance curves 
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FIG. 9: (color online): (a) PCAR spectra obtained in Bai_xKxFe2As2 single crystals by injecting the current along the ab 
plane |56l 157) and along the c axis [59]. The curves are plotted as a function of eV/ksT^ because the contacts had different 
T^. Dashed lines represent the normal-state conductance at T^. (b) The same curves after normalization (symbols) and 
the relevant fit with a two-gap 2D BTK model (lines). The values of the gaps and of the weight wi are reported in labels. 
Vertical lines indicate the position of the gap structures. Inset: temperature dependence of the gaps measured along the ab 
plane [561 157) . (c) Low temperature curves measured in Ba(Fei_xCox)2As2 single crystals, with current along c (curves 1-4) 
and ab (curves 5 and 6). From |56j (curves 1-2) and )22| (curves 3-6). Thin (thick) lines represent the 1-gap (2-gap) 2D BTK 
fit. Arrows indicate bosonic structures, (d) Experimental conductance curves of a c-axis contact as a function of temperature, 
with the relevant 2-gap fit. The gap values are reported in the inset. 



with no Andreev signal when the tip was gently put 
in contact with the cleaved sample surface, and a 
zero-bias conductance enhancement when the pressure 
was increased sufficiently. This feature is clearly related 
to superconductivity since it disappears at Tc, but its 
shape cannot be reconciled with an isotropic gap. One 
could think that this feature is rather suggestive of a 
d-wave order parameter, and try to fit it with a suitable 
generalized model. However, the current is here injected 
along the c axis. In this configuration, no interference 
effects are expected for the d-wave symmetry (with nodes 
on the hole FS), for the so-called "nodal s" symmetry 
(with nodes on the electron FS) [3] and even for the 
particular 3D nodal symmetry specifically predicted for 
122 compounds More exotic symmetries (if any) 

should be imagined to make HLQ and ELQ experience a 
sign change of the pairing potential in this configuration. 
Alternatively, one could imagine the zero-bias maximum 
to simply arise from the anisotropy of the OP but in this 
shown in jll C[ a closed 3D FS has to be invoked, 
which is not predicted by bandstructure calculations. 

We recently performed PCAR measurements with the 
"soft" technique in single crystals of Ba(Feo.9Coo.i)2As2 
(bulk r°" = 24.5 K, STc = 1 K). We obtained various 
series of spectra by injecting the current along the basal 
plane or along the c axis [52]. The normal-state conduc- 
tance at presented a shallow V-shape and progres- 
sively flattened on increasing temperature. Figure |9jc) 
reports some examples of normalized (i.e. divided by the 
conductance at T^) and symmetrized curves for c-axis 



(curves 3,4) and a6-plane contacts (curves 5,6). They all 
feature maxima related to a small gap, plus additional 
shoulders that can be related to a second, larger gap. 
The position of these features does not depend on the 
direction of current injection and this speaks in favor 
of isotropic gaps. In this concern, it is worth mention- 
ing that also the relative weight of the two terms in the 
conductance (let's say, wi and (1 — wi)) is the same ir- 
respective of the direction of the current. This seems 
to indicate that all the FS sheets have a similar degree 
of three-dimensionality, as also shown by ARPES [ST] 
and x-ray Compton scattering |62| . Thick lines in Hg- 
ure^c) represent the two-gap 2D BTK fit of the curves, 
while thin lines indicate the single-gap 2D BTK fit. It 
is clear that the latter model is unapt to reproduce the 
shape of the curves. The two-gap model works better in 
the gap region but leaves out additional structures (indi- 
cated by arrows) which are very clear only in curves with 
a large amplitude and that can be related to the strong 
electron-boson coupling (see sections II D and IV). Inci- 



dentally, note that in curve 5 these structures are proba- 
bly masked by small dips, that however do not affect the 
clear double-gap structures. The absence of zero-bias 
peaks certainly rules out, also in this case, the presence 
of nodes or line zeros on the Fermi surface, and this is 
consistent with theoretical predictions [5j. Instead, these 
findings are not incompatible with the presence of gap 
minima or even zeros (but not line zeros) in some regions 
of the Fermi surface, as recently proposed to explain the 
Raman scattering data in this compound [63) . 
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The gaps obtained by PCAR can now be compared 
with the results of ARPES measurements in the super- 
conducting state. Terashima et al. [M] measured a large 
isotropic gap A2 ~ 7 ± 1 meV on the hole-like FS sheet 
around F and a smaller isotropic gap A2 ~ 5 ± 1 meV on 
the hybridized electron-like FS sheets around M. The val- 
ues of the gaps as a function of temperature are shown as 
open circles in the inset to figure [9]jd) . As for the values 
of Al, ARPES and PCAR agree rather well, while the 
value of A2 is rather different, although the error bars 
overlap. A possible explanation of this difference can be 
found in [55]. Also note that the ARPES gaps show a 
weaker temperature dependence and go to zero abruptly; 
the same result has been found in Bai_xKxFe2As2 [54] . 



3. Ca(Fei^^Co^)2As2 

Very recently, we have performed PCAR measure- 
ments on Ca(Fei_2,Coa;)2As2 single crystals, with a bulk 
critical temperature Tc — 18 K as determined from DC 
susceptibility measurements. The "soft" point contacts 
were made on the freshly exposed side surface of a single 
crystal. Figure 10 a) shows the temperature dependence. 



from 2.21 K up to 16.16 K, of the conductance curves of 
a a6-plane point contact (Rn — 3 f2) normalized to the 
normal-state conductance at — 18 K. The conduc- 
tance curves are very different from those observed in the 
other compounds presented so far. First of all, they have 
very small amplitude; then, they systematically present 
(in 100 % of the cases) a zero-bias enhancement of the 
conductance, similar to that occasionally observed for ex- 
ample in Bai_xKxFe2As2 [59] and in Nd-1111 [34]; third, 
they have a rather broad, almost triangular shape. These 
characteristics do not change sensibly when a magnetic 
field up to 6 T is applied parallel to the ab plane, as shown 
in panel (b). This rules out magnetic scattering and in- 
trinsic Josephson junctions as the origin of the peak and 
indicates that the upper critical field is much higher than 
the maximum applied field. Furthermore, the systematic 
occurrence of the ZBCP, the similarity between different 
series of curves obtained in contacts with different re- 
sistance and the absence of effects normally observed in 
non-ballistic contacts, point towards an intrinsic origin 
of the ZBCP. 

As shown in section [II C| a zero-bias peak or maximum 
in a6-plane contacts can be a sign of a nodal OP. Let's 
assume that this is the case. This could be in agreement 
with the predictions by Suzuki et al. [F that, in 122, 
nodes appear on the quasi-3D hole-like FS around the F 
point (called ai in [5 ) when the height of the pnicto- 
gen on the Fe layer is small. In the similar compound 
Ba(Fe, Co)2As2 [5S] bandstructure calculations show in- 
deed that one of the holelike cylinders is strongly warped 
and that its diameter grows considerably going from F 
to Z, which is the situation envisaged in [5]. How- 
ever, in the absence of specific data about the FS of 
Ca(Fei_2;Coa;)2As2 and of detailed information on the 



fine structural details, we must rely on the PCAR spec- 
tra alone. In the absence of a specific way to express the 
peculiar nodal symmetry predicted for 122 compounds, 
we just tried to see whether a simple nodal symmetry 
(i.e. the d-wave one) can account for the observed PCAR 
results. Note that quite similar results can be however 
obtained by fitting with a fully anisotropic s-wave OP as 
the one discussed in section [Tl C[ To reproduce the ZBCP 
a misorientation angle a between the lobes of the d-wave 
gap and the normal to the interface has to be included 
in the 2D BTK model (which is the simplest approxi- 
mation, justified by the roughness of the fit). We were 
able to fit the conductance curves at any temperature 
up to = 18 K by using a small c?-wave gap and a 
larger s-wave gap (presumably associated to the electron 
FS sheets). As usual, we kept the barrier parameters 
Zi — 0.5 and Z2 = 0.285 as well as the weight wi — 0.55 
and the angle a = n/Q independent of temperature. The 
fit of a subset of the curves in figure 10 [a) is shown in fig- 
ure 10 c), and the relevant gaps are shown as a function 



of temperature in panel (d). Although the statistics is 
still to be completed, these results may suggest the Co- 
doped Ca-122 system as a possible good candidate for the 
study of nodal (or fully anisotropic) superconductivity in 
122 compounds. 



C. "22426" compounds 

Evidences of nodal OP have also been reported for 
another class of compounds, the highly layered "22426" 
family. Yates et al. [55] performed PCAR experiments, 
either with the conventional "needle-anvil" technique or 
the "soft" technique in Sr2ScFeP03 polycrystals with 
Tc ~ 17 K. This material is particulary interesting, 
though less known than other Fe-based compounds, for 
the study of the OP symmetry. As a matter of fact, in 
fill and 122 fS] the height of the pnictogen (here P) 
above the Fe plane has been indicated as the responsible 
for the appearance of nodes in the OP (in fill com- 
pounds, this is accompanied by a decrease of Tc which 
does not necessarily occur in 122). In Sr2ScFeP03 the 
angle of the Fe-P-Fe bond (118°) and the a-axis lat- 
tice parameter (4.016 A) mean that the height of the 
P above the Fe plane is 1.20 A, which lies between that 
of LaFeAs(0,F) and that of LaFePO (in which some ev- 
idences of nodal superconductivity have been found). 

The spectra obtained in [6^ systematically show a 
zero-bias peak irrespective of the junction resistance 
(ranging from 30 to 100 fl) and of the PCAR technique 
'a) shows the temperature dependence of 
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used. Fig 

one of the spectra obtained with a tip pressed against 
the sample surface. The curves are not offset; the ap- 
parent shift is due to a large conductance excess that 
extends up to very high energies. The inset shows indeed 
that it is only at eV > 100 meV that the tails of all the 
curves are superimposed. It is also clearly seen that the 
normal state conductance at T^ is approximately flat. 
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FIG. 10: (color online): (a) PCAR spectra obtained in 
Ca(Fei_xCox)2As2 single crystals by injecting the current 
along the ab plane, normalized to their normal state at 
— 18 K. (b) Magnetic-field dependence of the same spec- 
tra, at T = 2.2 K. The inset shows some raw conductance 
curves at different temperatures, (c) Some of the spectra of 
panel (a) with the relevant two-band 2D BTK fit with a large 
isotropic gap A2 and a small d-wave gap Ai. (d) Tempera- 
ture dependence of the gaps as obtained from the fit of panel 
(c). Lines are BCS-like curves. 



If this curve is used to normalize the low-temperature 
conductances, the resulting curves will look very simi- 
lar to the raw ones, apart from the vertical scale. It is 
clear that a fit of these curves in the whole voltage range 
is impossible with any version of the BTK model with 
energy-independent order parameters. 

Figure [TI|b) shows the magnetic-field dependence of 
the conductance curves (unfortunately, of a different con- 
tact) at T = 9 K, normalized at 20 meV The 
ZBCP is progressively depressed by the field, as well as 
the overall Andreev signal. A relatively small field of 7 
T is sufficient to completely suppress superconductivity. 
Sr2ScFeP03 is thus peculiar since the normal state con- 
ductance at T < is experimentally accessible and can 
in principle be used to normalize the data. The inset re- 
ports indeed the raw conductance curve at 9 K in zero 
field (top) and at the upper critical field (bottom). The 
latter is different from the normal state at shown in 
the inset to panel (a) but unfortunately the two sets of 
curves refer to different contacts and it is not possible to 
draw conclusions in this regard. 

A detailed analysis of the field-dependence |66| shows 
that the ZBCP is likely to have an intrinsic origin, apart 
from a small contribution from intergrain weak links 
which is suppressed already at 1 T. However, the at- 
tempts to fit the conductance curves with a single d-wave 
gap were not completely satisfactory due to the residual 
excess conductance at high bias (that persists even if one 
uses the normal-state conductance at Hc2 to normalize 
the curves). Even a two-gap model fails unless abnor- 
mally high values of the large gap are used to simulate the 
excess conductance. Only when choosing a normalization 
that removes this excess, the authors obtained a fairly 
good fit of the central peak that gave A = 4.34 ± 0.04 
meV, corresponding to a gap ratio 2A/kBTc ~ 6.7. 
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FIG. 11: (color online): (a) PCAR spectra obtained in a 
Sr2ScFeP03 polycrystal by Yates et al. |66] at different 
temperatures. The curves are not offset. Inset: the low- 
temperature raw conductance (top) and the normal-state one 
at (bottom), (b) Magnetic-field dependence (at 9 K) of 
the conductance curves of a different point contact. All the 
curves were normalized to their value at 20 meV. Inset: the 
raw conductance in zero field (top) and at Hc2 (bottom). 



IV. ELECTRON-BOSON INTERACTION IN 
FE-BASED SUPERCONDUCTORS 

As explained in detail in section |II D[ in moderate- 
or strong-coupling superconductors the PCAR conduc- 
tance curves can show signatures of the energy depen- 
dence of the superconducting gap. These structures are 
more easily observable if the amplitude of the Andreev 
signal is large; their signature in the second derivative of 
the I — V curve can be related to the electron-boson spec- 
tral function. In the following, we will show some results 
obtained in Ba-122 single crystals and Sm-1111 polycrys- 
tals that provide examples of such strong-coupling effects 
and of their analysis. Especially for the case of Ba-122 
compounds, when compared to inelastic neutron scatter- 
ing measurements, these results strongly support a spin- 
fiuctuation-mediated origin of superconductivity. 
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A. "122" compounds 

Figure [T2{ ^ a) shows an experimental conductance curve 
(circles) obtained on a Co-doped Ba-122 single crystal 
close to optimal doping. It can be clearly seen that the 
curve features additional structures at energies higher 
than the large gap, around 20 mV and, (less clearly) 
at about 40 mV. The solid line represents a theoretical 
conductance curve obtained by introducing in the two- 
band 2D BTK model the energy-dependent supercon- 
ducting gaps calculated within a three-band Eliashberg 
model (further details can be found in [55j) by using a 
Lorentzian electron-boson spectrum with characteristic 
boson energy f^o = 12 meV, very similar to the spin res- 
onance energy observed in neutron scattering measure- 
ments j31j . The theory can reproduce the shape and po- 
sition of the additional structure at 20 mV while the one 
at higher energy is only faintly visible in the theoretical 
conductance curve. Figure 12 b) shows the temperature 



dependence of the sign-changed second derivatives (solid 
lines) of the I-V curves of the same contact. As it can 
be expected, the electron-boson interaction features are 
now seen more clearly. The dashed line is the derivative 
of the theoretical curve shown in figure 12 a): it repro- 
duces very well the 20 mV-structure while the other one, 
though present, is hardly visible. This could be due to 
the fact that the actual shape of the electron-SF spec- 
tral function is more complex than the simple Lorentzian 
used here. By looking at the temperature dependence of 
the —(Pl/dV^ curves (another example is reported in fig- 
ure 12 ^c)) it is possible to notice that, as expected, both 
structures shift in energy on increasing temperature. The 
position of the lower-energy peak as a function of tem- 
perature, Ep(T) is reported in figure [T2|^d) for the two 
cases (full symbols) . The decrease of Ej, is at least par- 
tially due to the closing of the superconducting gap with 
increasing temperature. It is thus essential to extract 
the temperature dependence of the characteristic energy 



of the boson spectrum, fltiT) 
section 



Ep{T) - A 

7nax 

(T) (see 



II D). Figure [T2|^d) shows that the values of Qb 
(open symbols) decrease on increasing temperature, and 
tend to zero when T T^. This indicates that Q(, cannot 
be the energy of a phonon mode (in that case it would 
not tend to zero!) and thus rules out a phononic origin 
of this feature. Instead, the trend of fij, is very similar to 
that of the spin-resonance energy peak reported in [ 31] 
and thus strongly supports a spin-fluctuation-mediated 
pairing mechanism in these compounds. 

In this regard it is worth noticing that the existence of 
a relationship between Tc and the spin-resonance energy 
observed by inelastic neutron scattering in Fe-based com- 
pounds does not necessarily prove the spin-fluctuation 
pairing. As J. Paglione and R.L. Greene point out, since 
such a relationship appears exclusively in nearly mag- 
netic unconventional superconductors, it could be consid- 
ered an independent remnant of nearby magnetic order 
rather than a strong signature of magnetically mediated 
pairing |67) . 



What PCAR spectroscopy provides is something more 
than that. Indeed: i) it detects features related to the in- 
teraction of electrons with elementary excitations in the 
material (and not only to the excitation itself); ii) it 
shows that these excitations enter the imaginary part of 
the energy-dependent order parameter; iii) it allows ex- 
tracting the characteristic energy of the interaction spec- 
trum a^F(Jl), which turns out to coincide with that as- 
sociated to spin fluctuations and observed by neutron 
scattering experiments. 



"S 1.3 



1.2 



!1.0 



(5 



> 
■a 



O Exp 
Eiashberg m \ 


(a) 


-40 -20 


20 40 


Voltage 


(mV) 


:(o) 


'L4.3 K 


Y \ ^ 8.0 K 




,^ 12.2 K 

15.2 K 




17.1 K 


19-3 K 




10 20 30 40 


Voltage 


(mV) 




10 20 30 

Voltage (mV) 



40 




Temperature (K) 



FIG. 12: (color online) (a): Normalized experimental conduc- 
tance curve (circles) obtained in a Ag/BaFei.8Coo.2As2 point- 
contact together with the theoretical one (line) obtained by 
introducing in the BTK model the energy-dependent super- 
conducting gaps. The gaps as obtained by the Eliashberg cal- 
culations are Ai — 6.1 meV, A2 — —3.8 meV and A3 — —8.0 
meV [22]. (b): Temperature dependence of the —cfl/dV^ 
curves obtained from the same contact of (a) showing the 
displacement of the bosonic structures with increasing tem- 
perature. Dashed line is obtained from the theoretical curve 
shown in (a), (c): same as in (b) but for another contact, 
(d): Temperature dependence of the energy peak, Ep (full 
symbols) extracted from (b) and (c) together with the corre- 
sponding boson energy Qb{T) = Ep — /S.max- Lines are guide 
to the eye. 



B. "1111" compounds 

Figure [Tsj^a) shows a normalized experimental conduc- 
tance curve obtained by us on a optimally doped Sm-1111 
polycrystal (circles). It looks very similar to the spec- 
trum by Naidyuk et al [JT] shown in figure [sj^d) and also 
to some PCAR spectra in underdoped polycrystals by 
Chen et al. [52]. The amplitude of the Andreev signal 
in this contact is exceptionally high (about 80%) and, 
besides clear two-gap features (peaks and shoulders) , ad- 
ditional structures or small kinks can be seen around 27 
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mV and 40 mV. The dashed hne is a BTK fit to the exper- 
iment by using the two-band 2D BTK model with ECS 
gap values (i.e. independent of energy). As in figurejsjd), 
the fit reproduces very well the experiment in the central 
part of the curve (and allows obtaining reliable values of 
the gaps) but fails at higher energies. The solid line is in- 
stead the result of inserting in the same BTK model the 
energy-dependent OPs obtained by solving a three band 
Eliashberg model [3S] in which, as usual, the electron- 
boson spectral function is modeled by a Lorentzian curve 
and the calculated gaps are Ai = 17.23 meV, A^6.03 
meV and A3 — —19.56 meV. 



■ r^o), the feature at 40 mV is remarkably well re- 



produced as can be observed in figure [T3[a) and 13 b) 
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FIG. 13: (color online) (a): Normalized experimental conduc- 
tance curve (circles) obtained in a Ag/SmFeAsOo.8Fo.2 point- 
contact. Dashed line is a BTK fit to the experiment, obtained 
by using the constant BCS values for the gaps. The param- 
eters of the fit are Ai = 6.0 meV, A2 = 19.5 meV, Fi = 0.1 
meV, r2 = 2.5 meV, Zi=0.2, Z2=0.1 and wi=0.6. Solid line 
is a theoretical curve obtained by introducing in the BTK 
model the energy-dependent gap functions calculated within 
the three-band Eliashberg theory whose values are Ai = 6.03 
meV, A= 17.23 meV and A3 = -19.56 meV. (b): Tempera- 
ture dependence of the —d?I/dV^ curves obtained from the 
same contact of (a) showing the displacement of the bosonic 
structures with increasing temperature. Dashed line is ob- 
tained from the theoretical curve shown in (a). Inset: Tem- 
perature dependence of the energy peak, Ep (full symbols) 
extracted from (b) together with the corresponding boson en- 
ergy Q.b(T) = Ep{T) — A,nax{T). Lines are guide to the eye. 

Since to the best of our knowledge no spin-resonance 
energy value is available for this compound, the char- 
acteristic energy has been chosen according to JQJ\ by 
extrapolating the relationship flo ^ A.GdksTc ^ 20 meV. 
Although the theoretical curve shows no structures at 
27 mV (which by the way corresponds approximately to 



(dashed line). Similarly to the case of Co-doped Ba-122, 
only the structure present at approximately f2o + ^max 
is reproduced. This again indicates that the model has 
to be investigated further or that additional features of 
the spectral function are playing an important role. As 
expected, both structures shift in energy on increasing 
temperature, partly because the amplitude of the su- 
perconducting gaps is also decreasing. The inset to fig- 
ure |13Ka) reports the position of the energy peak in the 



second derivative, Ep (full symbols) and the values of 
fib = Ep — Amax (open symbols) as a function of tem- 
perature. Similarly to the case of Co-doped Ba-122 re- 
ported in the previous section, Qi, decreases in tempera- 
ture indicating that the observed structure does not have 
a phononic origin and might instead be related to spin 
fluctuations. It is interesting to note that similar struc- 
tures at about 20, 40 and 55 meV were directly observed 
in the PCAR spectra by Chen et ai, who also observed 
their shift in energy, similar to that of the gap [52], but 
did not give them any special physical meaning. 



V. THE PUZZLE OF 2A/kBTc RATIOS 

Figure [14] reports a collection of the values of the gap 
ratios 2Ai/kBTc (where A^ is the i-th band gap measured 
at low temperature) given by different PCAR measure- 
ments performed by us and other groups in various Fe- 
based systems. Only results referring to nodeless gaps 
have been reported essentially i) because they represent 
the great majority the cases studied up to now and ii) 
to allow an homogeneous comparison with the standard 
s-wave BCS ratio. The gap ratios are reported as a func- 
tion of the Tc of the sample or, when available, of the 
of the contact. Two representative ARPES data are 
also shown for comparison. Let us start by analyzing the 
trend of 2Ai/kBTc, where Ai is, as usual, the small gap. 
In samples with high Tc, this ratio is close to (but gen- 
erally smaller than) the BCS value for the s-wave weak- 
coupling regime. In most compounds, it is around 2.5. 
It then remains almost constant when the Tc decreases 
down to approximately 30 K; then, it starts to increase 
continuously with further decreasing Tc- 

A qualitatively similar trend is also observed for 
2A2 /ksTc, despite a larger spread of values which is not 
surprising since the larger gap is worst resolved than the 
small one, its absolute value could be affected by the 
choice of the normalization and moreover, as pointed out 
in |IV[ could also be the average of two large gaps of sim- 
ilar value. We cannot, however, exclude that the spread 
is also due to specific features of each system such as 
nesting, vicinity to magnetism etc. but also, and most 
probably, quality of the sample, homogeneity and disor- 
der. 

Although more data points are desiderable in the low- 
Tc region of the graph, figure 14 indicates that the in- 
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FIG. 14: (color online): 2Ai/kBTc and 2A2/kBTc values vs 
Tc as determined by PCAR spectroscopy measurements per- 
formed by our group as well as several other groups in case of 
nodeless gaps. Horizontal dashed line represents the s-wave 
weak coupling limit for the BCS theory. Blue and red dashed 
lines are only guides to the eye. PCAR data for SmOFe(As,P) 
are taken from [68| . other results published in literature are 
from [33H361 BIl ET] El |66l 113 [ZQ] . The data labelled "Wang 
et al." indicate the two-band s-wave fit of the spectra from 
ref. [7T]- ARPES data are from E1|M]- Data from flTj are 
plotted vs T/. 



crease of the 2Ai/kBTc ratio below a certain Tc (on the 
order of 30 K) appears to be a rather common trend for 
both gaps. But, how can this trend be interpreted? 

In order to make a reliable statement in this regard, 
the complexity of these systems has to be carefully 
taken into account. Because of the multiband nature of 
these compounds and their dominant interband pairing, 
a proper interpretation of the data would require a quan- 
titative analysis (for example within a multiband Eliash- 
berg model) taking into account the electronic structure 
of each compound. However, as a first approximation, 
we have done numerical simulations within the three- 
band s-± Eliashberg theory (with feedback effect [28] . 
i.e. by including in a self-consistent way the temper- 
ature dependence of the spin fluctuation spectral func- 
tion) to see in what conditions this trend can be re- 
produced. To simplify the problem, we fixed some pa- 
rameters to the values that allowed us to reproduce the 
gaps in Sm-1111 compound with Tc = 52 K [55|. These 
parameters are A23/A13 (where Xij is the coupling con- 
stant between bands i and j), ni3 and ^23 {riij = Ni/Nj 
where Ni is the density of states of band i). Accord- 
ing to |67| . we assumed D,q to vary linearly with Tc, i.e. 
Hq ~ A.GBksTc- For values of the critical temperature 
down to 30 K, the experimental values of 2Ai/kBTc can 
be reproduced by keeping the total coupling constant, 
Xtot, practically constant to a value of approximately 2.4. 
Below 30 K, to follow the experimental trend it is nec- 
essary to i) increase Xtot and ii) modify the temperature 
dependence of the peak energy. In standard situations, 
the latter evolves with temperature according to fio (T) = 



rio tanh[1.76v'(Tc* /T - 1)] [2H] (i-e. in a BCS-like way) 
where T* is the critical temperature calculated without 
the correction due to the feedback effect (T* > Tc). Here, 
we assume instead n„{T) = fla tanh[1.76fc^(Tc*/T - 1)] 
where fc = 1 for Tc > 30 K and decreases for Tc < 30 K. 
Roughly speaking, this is equivalent to a reduction of the 
representative frequency Qq at any T between and T* , 
with respect to the above BCS-like curve. As a result, a 
larger Atot is necessary to obtain the correct Tc; but this 
also gives larger gaps. This result, if confirmed by using 
detailed electronic calculations for each compound, might 
indicate that below Tc=30 K the coupling increases and 
would also have interesting implications on the tempera- 
ture dependence of the spin fluctuation spectral function 



that, as shown in flgure 12 ^d) and in the inset of flgure 13 
can be determined from high-quality PCAR spectroscopy 
experiments. 



VI. SUMMARY AND CONCLUSIONS 

In the present Report we have presented an overview of 
the (partly) solved and of the (still) debated points about 
superconductivity in Fe-based superconductors that have 
been (or can be) addressed by means of point-contact 
spectroscopy. To do so, we have also introduced general- 
izations of the usual models for the fit of the PCAR spec- 
tra, and presented some new results. Here, we would like 
to simply summarize the main conclusions of this long 
discussion, and outline the possible future challenges for 
PCAR spectroscopy in the research about Fe-based su- 
perconductors. 

In Section |ll] we have examined and discussed in de- 
tail the profound and often counterintuitive connections 
among the FS topology of multiband superconductors, 
the symmetry of the CPs present on the FS sheets, the 
directionality of the PCAR process, the strength and the 
spectrum of electron-boson interaction, and the resulting 
PCAR normalized conductance curves. Taking into ac- 
count these connections greatly complicates the descrip- 
tion of the AR process, but is unavoidable when deal- 
ing with Fe-based superconductors. The extension of the 
BTK model to its most general form, that we call "full 
3D BTK" model, (sections [!lC| and |IID| ) allows directly 
comparing the PCAR experimental results with the pre- 
dictions of the theory for a specific shape of the FS, a 
specific symmetry of the OP in every band and a specific 
spectrum of electron-boson interaction (via the prelimi- 
nary solution of the Eliashberg equations). The compar- 
ison of the results to those of the simplified BTK models 
provides some indications and warnings for the analysis 
of PCAR data in Fe-based compounds that can be sum- 
marized as follows: 

1) The analysis of the a6-plane PCAR spectra on Fe- 
based superconductors requires at least the 2D version of 
the BTK model, because of the quasi-2D shape of their 
FS and the possible presence of non-isotropic OPs. Even 
in the case of pure s-wave symmetry of all the OPs, this 
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model gives results slightly different from the full 3D one 
(in particular, an apparent enhancement of the Z values 
is observed, due to FS-shape-related geometric effects); 

2) In principle, c-axis PCAR spectra in pnictides 
should be analyzed only by means of the 3D BTK model 
(simplified version of [TS] or full version of equation 11) 
because of the non-spherical shape of the FSs. However 
if the a&-plane contacts on the same material don't show 
any sign of line nodes or zeros in the OP, and the precise 
determination of the Z values is not vital, the 2D BTK 
model can return correct information about the isotropic 
gap values and their temperature dependency. In this 
last case, as shown for the large gap in figure |4]^b), the 
Z-enhancement effect is remarkable especially when the 
FS is only slightly warped, but the gap-related peaks re- 
main at the same energy positions; 

3) The apphcation of the full 3D BTK model to the 
case of fully anisotropic s OPs (with line zeros) or rf-wave 
OPs (with line nodes) that open on warped cylindrical 
surfaces leads to interesting conclusions. As expected, 
the normalized conductance of a&-plane contacts in the 
AR regime {Z < 0.2 — 0.4) shows characteristic zero- 
bias maxima or peaks, respectively. Unfortunately, in 
the presence of a small lifetime broadening the smoothed 
spectra become very similar, preventing the possibility to 
discriminate between the two symmetries. In the same 
conditions, c-axis contacts unexpectedly show no signs 
of zero-bias peaks or maxima even for Z = because of 
the geometric Z-enhancement effect, as shown in figure [4] 
This prediction could be verified by means of directional 
PCAR experiments on high-quality single crystals. In 
other words, the observation of zero-bias peaks or max- 
ima in c-axis contacts can be reconciled with the previous 
two OP symmetries only if the corresponding FS sheet 
is 3D, thus giving to directional PCAR spectroscopy the 
capability to attain indirect information on the shape of 
the FS; 

4) The introduction of energy-dependent gap functions 
Ai{E) (determined by the proper solution of multiband 
Eliashberg equations) into the BTK expressions leads to 
the conclusion that structures due to the electron-boson 
interaction (EBI) are certainly visible in the AR con- 
ductance at energies higher than the large gap (and for 
coupling constants compatible with the Tc and the gap 
values) provided that the spectral function of this interac- 
tion is sufficiently peaked (e.g. has a Lorentzian shape). 
The use of the spectral function typical of antiferromag- 
netic spin fluctuations (AFSF) results in curves with no 
EBI structures unless unphysical, very low cutoff energies 
are chosen. Serious questions thus arise on the true shape 
of the bosonic spectrum since (as shown in two particu- 
lar cases in section IV I EBI structures have actually been 



observed in the PCAR data of various compounds. 

As for the OP symmetry in Fe-based superconductors, 
the analysis and comparison of various PCAR results re- 
ported in section [lTT| show that, despite the rather com- 
mon opinion that PCAR results are contradictory, the 
most recent and reliable experiments are in fairly good 



agreement with one another. In particular: 

1) When data from various groups are available for a 
given compound, they generally give similar results if an- 
alyzed in a homogeneous way. The apparent conflicts are 
mostly due to a different normalization of the raw data 
or to a different interpretation of the common structures 
present in the conductance. 

2) In optimally doped and slightly underdoped 1111 
compounds, if the PCAR curves are normalized to 
the relevant normal-state spectrum and the high-energy 
structures are not artificially reduced, all the data point 
to the presence of (at least) two isotropic gaps (with no 
nodes and no line zeros) with ratios 2Ai/kBTc lower and 
sensibly higher than the BCS value, respectively. 

3) Similar observations apply to the 122 compounds 
where data of various groups concur in supporting a 
multi-gap s± picture (though at most two gaps can be 
resolved by PCAR spectroscopy). Since in this case most 
of the experiments are made on single crystals, a depen- 
dency of the "weight" of every band on the direction of 
current injection is also observed, in agreement with the 



predictions of subsection II C 



4) In compounds of the 22426 family and in Co-doped 
Ca-122 (as shown here for the first time) PCAR spec- 
troscopy gives some evidence of a d-wave or a fully 
anisotropic s-wave OP. In these cases, experiments along 
different crystallographic directions can be used as a 
"test" of the predictions of II C or of the presence of a 
3D FS sheet. 



In section IV we have given evidence of EBI struc- 
tures in the PCAR spectra of Co-doped Ba-122 and Sm- 
1111. By properly combining the direct solution of the 
three-band s± Eliashberg equations and the BTK mod- 
els we were able to reproduce these structures in ampli- 
tude and position, thus obtaining important information 
on the characteristic energy of the bosonic spectrum and 
its temperature dependency. This energy turns out to 
coincide with the characteristic energy of the spectrum 
of spin fluctuations measured by inelastic neutron scat- 
tering thus giving a strong support to the SF origin of 
the coupling in these compounds. In addition these EBI- 
related structures (not only small peaks or shoulders at 
the proper energy but also a general increase of the con- 
ductance at high bias) can nicely explain some frequently 
observed anomalies of PCAR spectra that cannot be re- 
produced by the BTK models using constant gaps. This 
is another example that too simplified models cannot ex- 
plain the rather complex physics of AR in Fe-based su- 
perconductors. 

Finally, in section[V]we have shown ~ by comparing the 
results of various experiments of different groups - that at 
the decrease of the Tc of the sample the ratios 2Ai/kBTc 
for the small and the large gap have an anomalous be- 
havior. Instead of merging at a value close to BCS one at 
the lowering of Tc, as naively expected in the framework 
of a multiband model dominated by interband couplings, 
they both increase at T < 30 K reaching values greater 
than 5 and 10, respectively. A preliminary analysis of 
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this "strange" behavior in the framework of the Ehash- 
berg theory including feedback effect [28j shows that it 
can be explained only supposing that below ^ 30 K the 
temperature dependency of the characteristic boson en- 
ergy is no more BCS-like but becomes more linear and, 
at the same time, the total coupling strength increases. 
Both these conclusions could be further investigated by 
directional PCAR experiments in low- Tc high-quality sin- 
gle crystals. 

In conclusion, we have shown that PCAR spectroscopy 
can provide information on some debated properties of 
the Fe-based superconductors if the PCAR spectra are 
analyzed within suitably generalized models. In partic- 
ular, a tighter relationship between the electronic prop- 
erties (and thus the geometry of the Fermi surface) and 
the PCAR results can be established. This could help 
either verifying the plausibility of calculated Fermi sur- 
faces, or providing some experimental hints about what 
the FS should look like to explain the measured spectra. 
A systematic PCAR study of some selected systems as 
a function of doping, coupled to bandstructure and FS 
calculations as well as with the determination of fine lat- 
tice details, could help testing the proposed connection 
between these aspects and the symmetry of the order pa- 
rameter, and to verify to what extent the scaling law of 
the gap ratios l^ilksTc as a function of is univer- 
sal in these compounds (and, if so, what are its physical 



reasons and implications) . As for the origin of the cou- 
pling mechanism, important hints can be obtained by 
extracting the typical boson energy from PCAR spec- 
tra and comparing it with the energy of spin resonance 
as measured, for instance, by neutron scattering exper- 
iment. In general, thus, the future developments of the 
research about Fe-based compounds pose a challenge to 
PCAR spectroscopy, and might strongly benefit from the 
results of this experimental technique. 
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